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W)' Abstract 



We calculate the string tensions of fc-strings in SU(A^) gauge theories in both 3 and 4 
dimensions. We do so for SU(4) and SU(5) in D=3+l, and for SU(4) and SU(6) in 
D=2+l. In D=3+l, we find that the ratio of the k = 2 string tension to the k = 1 

CN I fundamental string tension is consistent, at the 2a level, with both the M(-theory)QCD- 

inspired conjecture that at oc sin(7rA;/A^) and with 'Casimir scaling', ak oc k{N — k). In 

O ■ D=2+l, where our results are very precise, we see a definite deviation from the MQCD 

^ ! formula, as well as a much smaller but still significant deviation from Casimir scaling. 

Q I We find that in both D=2+l and D=3+l the high temperature spatial fc-string ten- 

sions also satisfy approximate Casimir scaling. We point out that approximate Casimir 
scaling arises naturally if the cross-section of the flux tube is nearly independent of the 
c^ ■ flux carried, and that this will occur in an effective dual superconducting description, 

if we are in the deep-London limit. We estimate, numerically, the intrinsic width of 
D I /c-strings in D=2+l and indeed find little variation with k. In addition to the stable 

^. [ /c-strings we investigate some of the unstable strings, which show up as resonant states 

in the string mass spectrum. While in D=3+l our results are not accurate enough to 

^ ■ extract the string tensions of unstable strings, our more precise calculations in D=2+l 

show that there the ratios between the tensions of unstable strings and the tension of 
the fundamental string are in reasonably good agreement with (approximate) Casimir 
scaling. We also investigate the basic assumption that confining flux tubes are de- 
scribed by an effective string theory at large distances, and attempt to determine the 
corresponding universality class. We estimate the coefficient of the universal Liischer 
correction from periodic strings that are longer than 1 fermi, and find Cl = 0.98(4) in 
the D=3-|-l SU(2) gauge theory and cl = 0.558(19) in D=2+l. These values are within 
2(j of the simple bosonic string values, cl = vr/3 and c/, = 7r/6 respectively, and are in- 
consistent with other simple effective string theories such as fermionic, supersymmetric 
or Neveu-Schwartz. 

PACS Numbers: 11.15.-q, 12.38.Aw, 11.15.Ha, 12.39.Pn. 

Key Words: SU(iV) Gauge Theories, Lattice Gauge Theories, Confinement, String 

Tension, Strings in Higher Representations. 



1 Introduction 

It is widely believed that the SU(3) gauge theory that underlies QCD is linearly confin- 
ing and that this explains why we do not observe quarks (or gluons) in nature. The fact 
that confinement is linear suggests that the colour-electric fiux between fundamental 
charges is localised in a tube between those charges and it is attractive to think that 
the long-distance physics of such fiux tubes is given by an effective string theory. The 
simplest possibility is that this string theory is bosonic but other possibilities are not 
excluded and indeed might be natural if QCD is obtained by some kind of reduction 
from a higher-dimensional theory. 

The same comments apply to SU(A^) gauge theories for A^ 7^ 3. Indeed there are 
long-standing ideas that for iV — > 00 the SU(A^) gauge theory can be thought of as a 
string theory. Moreover SU(iV) gauge theories in D=2+l also appear to be linearly 
confining [|I| and all the above comments will apply there as well. 

In addition to charges in the fundamental representation (like quarks) one can con- 
sider the potential between static charges in higher representations of the gauge group. 
In SU(2) and SU(3) any such charge can be screened by gluons either to the funda- 
mental or to the trivial representation. Since virtual gluons are always present in the 
vacuum this means that such a potential will, at large distances, either rise linearly 
with a string tension equal to the fundamental one or will fiatten off to some constant 
value. (This assumes that the fundamental string tension, a, is the lowest, as appears 
to be the case.) For A^ > 4, however, this is no longer the case and there are new stable 
strings with string tensions different from the fundamental one. A typical source may 
be thought of as k fundamental charges located at a point. The confining string is then 
usually referred to as a /c-string. For SU(A^) we have non-trivial stable /c-strings up to 
a maximum value of k given by the integer part of N/2. 

Such /c-strings are interesting for a variety of reasons. The values of their string 
tensions, Ufc, will constrain models of confinement. In models of glueballs in which 
the latter consist of open or closed strings, the SU(A^) mass spectrum should change 
with A^ in a way that is determined by how ak varies with A^ and k. In addition there 
are theoretical ideas concerning the value of ak- In particular there is a conjecture 
based on M-theory approaches to QCD (MQCD) that suggests a/, oc sm{nk/N}. 
One can contrast this with the old "Casimir scaling" conjecture that would suggest 
CTjfc oc k{N — k) and also with the simple possibility that a fc-string consists of k non- 
interacting fundamental strings, in which case at = ka. 

In a string theory the mass of a fiux tube of length / will receive a leading large-/ 
correction that is 0(1//). Such a slowly decreasing correction cannot be made negligible 
simply by making / ^ Ifm and so it will, in principle, limit the accuracy of our 
calculations of ak- Fortunately this leading string correction is known to be universal 
[Q in that its coefficient is determined entirely by the central charge of the effective 
string theory. The universality class is usually thought to be that of a simple Nambu- 
Goto bosonic string. There is however no strong direct (numerical) evidence for this 



belief that we are aware of. Such evidence would need to be obtained from strings that 
are longer than 1 fm and to achieve the required accuracy for such strings is a hard 
numerical problem. Where accurate values are quoted in the literature they typically 
involve fitting potentials down to shorter distances, where the fits are almost certainly 
dominated by the tail of the Coulomb term which has the same functional form as the 
string correction (and in practice a similar coefficient). We have therefore attempted to 
provide a usefully accurate calculation of this string correction in SU(2) gauge theory, in 
both D=3+l and D=2+l. Such a calculation also addresses the fundamental question 
of whether a confining fiux tube is in fact described by an effective string theory at 
large distances. 

The contents of this paper are as follows. In the next Section we describe how we 
calculate at from the mass of a flux loop that winds around the spatial torus. We 
contrast this method with one that uses explicit sources; in particular for 'strings' that 
can break. All this requires a classification of strings in all possible representations 
of the gauge group (the details of which appear in Section [A.l| of the Appendix). In 



the following Section we summarise the lattice aspects of our calculation; we can be 
brief since it is entirely standard. We then turn to the basic question of whether we 
really do have strings and, if so, which universality class they belong to. Confining 
ourselves to fiux tubes that are longer than Ifm we find that in both D=3+l and 
D=2+l SU(2) gauge theory the leading correction to the linear dependence of the 
string mass is consistent, within quite small errors, with what one expects from the 
simplest effective bosonic string theory and excludes the most obvious alternatives. We 
then turn to our calculation of fc-strings. We begin by briefly summarising some of the 
theoretical expectations: MQCD, Euclidean and Hamiltonian (lattice) strong coupling, 
Casimir scaling, the bag model and simple flux counting. We then turn to our D=3+l 
calculations of ak=2 in both SU(4) and SU(5) gauge theories and follow this with our 
(inevitably) much more accurate D=2+l calculations for SU(4) and SU(6). (In SU(6) 
we are able to address non-trivial A; = 3 strings.) In D=3+l we flnd consistency at the 
2cr level with both MQCD and Casimir scaling. In D=2+l the string tension ratios, 
while still close to the MQCD formula, are much closer to Casimir scaling. We point 
out that if the flux is homogeneous, then (approximate) Casimir scaling arises if the flux 
tube width is (approximately) independent of k. (And, more theoretically, that this 
will arise in the deep-London limit of a dual superconducting vacuum.) To test this idea 
we perform an explicit calculation of the intrinsic size of fc-strings in D=2+l. We flnd 
that the fc-string width is indeed largely independent of /c, albeit with some interesting 
if weak differences. We then point out that the same calculations can be reinterpreted 
as telling us that the spatial string tension in the high temperature deconflning phase 
satisfles approximate Casimir scaling. We complement this with an explicit D=3+l 
high-T calculation that demonstrates that in that case too the string tension ratio is 
close to Casimir scaling. We then attempt to see if there is any sign of other, unstable, 
strings, which should appear as excited states in the string mass spectrum. We flnd 



reasonably convincing evidence for such strings, satisfying approximate Casimir scaling, 
in our D=2+l calculations. We finish with a discussion of our results and some of their 
implications. 

A preliminary version of our D=3+l calculations has appeared in ||^. In the intro- 
duction to a recent companion paper on the mass spectrum and topological properties 
of D=3+l SU(A^) gauge theories we briefly summarised some of our results on a^ja. 
In particular we drew attention to the relevance of these results on /c-strings for the 
Casimir scaling hypothesis. We remark that all these calculations are intended as a 
first step to a much more complete and accurate calculation of the properties of SU(A^) 
gauge theories for all values of A^. 

2 Strings and string breaking 

Consider a static source in some representation 7?, of the gauge group in, say, 3+1 
dimensions. Suppose we have a conjugate source a distance r away. If r is small then 
the potential energy will be dominated by the Coulomb term 

V^ir) '=' ^^^^ii^ + . . . (1) 

r 

where as{r) is the usual running coupling and C-n is the quadratic Casimir of the 
representation TZ: 

Cn = TrnT'^T- (2) 

with the T" being the generators of the group. If the theory is linearly confining, and 
if we ignore the fact that the source may be screened by gluons, then at large r we 
expect the potential energy to be given by 

Vnir) = anr + .... 3 

24 r 

Here a-ji is the 'string tension' of the confining flux tube joining the sources and how its 
value varies with the representation TZ is an interesting physical question. If the long- 
distance physics of the confining flux tube is described by an effective string theory, then 
the 0(l/r) correction in eqn(^ is the Casimir energy of a string with fixed ends, and Cs 
is proportional to the central charge. This correction is universal 0], since it depends 
only upon the massless modes in the effective string theory and does not depend upon 
the detailed and complicated dynamics of the flux tube on scales comparable to its 
width. The central charge is given |^ by the number of massless bosonic and fermionic 
modes that propagate along the string. In practice it is usually assumed that c^ = 1, 
corresponding to the simplest possible (Nambu-Goto) bosonic string theory. However, 
these modes are not related to the fundamental degrees of freedom of our SU(A^) 



gauge theory in any transparent way and the presence of fermionic modes is certainly 
not excluded. For example, we have the following simple possibilities [|[: 

bosonic 

fermionic , . 

supersymmetric 

Neveu — Schwartz 

Whether a string description of the confining flux tube is in fact valid and, if so, what 
is its universality class, are fundamental questions which are still largely open. The 
examples in eqn(^) show that one needs to calculate Cs to better than, say, ±15% if 
one is to usefully resolve different possibilities. 

In reality the vacuum contains virtual gluons which can screen the static source, and 
this will complicate any attempt to calculate a-ji. When and how this happens will 
depend on the energetics of the system. Suppose the representation TZ can be screened 
to a different representation TZ' by a number of gluons. Such a screened source will 
acquire an extra mass of, say, AM. If the string tension corresponding to TZ' is smaller 
than a-ji, the screening is certain to become energetically favoured for sufficiently large 
r, since as r — > oo 

AM < Vn{r) - \/7^'(r)~(77^r - a^'r. (5) 

The minimum value of r at which the energetics favours screening is the string breaking 
scale Tft. (We shall use the term 'string breaking' when a source is screened to a different 
representation, even if the latter is not the trivial one.) So if we calculate the potential 
for our sources we can expect the r-dependence to be given by VTi{r) for r < Vf, and 
Vjz'ir) for r > Vf,. If r^ is large enough then we will be able to extract a-ji from the 
linearly rising potential at r < r^. In practice, however, the string breaking scale is 
similar to other dynamical scales in the theory and it is not clear whether any apparent 
linear rise of V{r) for r < r?, is due to the precocious formation of a string, from which 
we can read off a-ji, or if it is merely accidental. Indeed it may be that one cannot assign 
an unambiguous meaning to the quantity a-ji under these circumstances. However it 
is also possible that if the string breaking is relatively weak then one may be able to 
calculate a-ji for r > r^ by identifying an appropriate excited string state. In any case 
it is clear that string breaking creates substantial extra ambiguities in any attempt to 
calculate the properties of strings corresponding to higher representation charges. 

For SU(2) and SU(3) any representation TZ can be screened by gluons to either the 
trivial or the fundamental representation. However for SU(iV > 4) this is no longer 
the case and one finds new strings that are completely stable and to which none of 
the above ambiguities apply. The situation may be summarised as follows. (We leave 
a fuller discussion to the Appendix.) Suppose the representation 7Z can be obtained 
from the product of n-n + k fundamental representations and n-ji conjugate ones. Let 
z be an element of the centre, Z^, of the SU(A^) gauge group. Under such a centre 



gauge transformation the source will transform as z^. We shall refer to k as the J\f- 
ality of the representation. Now, since gluons transform trivially under the centre, 
the source will continue to transform in this way even if it is screened by gluons to 
some other representation. Thus the same value of k will label a source and all the 
sources that can be obtained from it by screening. Indeed one can show that if two 
representations have the same value of k then either one can be screened by gluons to 
the other. Within any given class of such sources there will be a lowest string tension 
ak, which, by string breaking, will provide the potential for any of these sources at 
large enough distances. The independent values of k are constrained: under charge 
conjugation k -^ —k, and we also have z^ = z^~^ . Thus for SU(A^) we have stable 
strings labelled by k = 1, ..., kmax where kmax is the integer part of N/2 and k = 1 is, 
of course, the fundamental string. That is to say, we must go to at least SU(4) to have 
a. k = 2 string, and to at least SU(6) to find a. k = 3 string. 

In this paper we shall compare the k = 2 and k = 1 string tensions in SU(4) and 
SU(5) gauge theories in D=3+l. We shall do the same in SU(4) and SU(6) in D=2+l; 
and in this last case we shall calculate the k = 3 string tension as well. In all cases 
we shall extrapolate to the continuum limit and the aim is to obtain results that are 
accurate enough to distinguish between various theoretical expectations. Since these 
strings are all stable there is no intrinsic ambiguity in defining a string tension and we 
can, in principle, achieve this goal. 

We shall calculate ak not from the potential between static charges but from the 
mass of a /c-string that winds once around the spatial torus. If the string length / is 
sufficiently large, its mass will be given by an expression similar to eqn(^): 

/^oo , n{D-2) cs 

mk[l) = (Tkl - + .... (6) 

o l 

We note that because of the different boundary conditions on the ends of the string 
(periodic rather than fixed) the 0(l/r) universal string correction is four times as large 
as for the static potential . We further note that because there are no explicit sources 
there is no analogue, at small /, of the Coulomb potential in eqn(|l]). That is to say, 
this is a particularly favourable context in which to calculate the string correction: its 
coefficient is large, and there is no danger of confusing it with a Coulomb interaction 
which has the same functional form. 

One can of course consider such closed but non-contractible winding strings in any 
representation 7^. However, just as with the static potential, such a string can be 
screened to a different string, corresponding to a representation TZ', as long as both 
strings possess the same A/"-ality. One can picture the string breaking as follows: a pair 
of gluons pops out of the vacuum somewhere along the string. These then move away 
from each other along the string. As they do so the section of string between them will 
no longer belong to TZ but rather to the product of 71 and the adjoint representation. 
If the two gluons propagate all the way around the torus they can meet and annihilate 
leaving a new string that is entirely in this different representation. Clearly one can 



extend this to any number of gluons. This is just hke the breaking of the string between 
static sources except that here the gluons eventually annihilate rather than adhering 
to a source. Thus there is no extra mass AM to consider and the breaking can occur 
for small / if this lowers the mass of the loop. That is to say, there is no region r < Vf, 
where one might hope to see a portion of the original string prior to its breaking. Of 
course, just as for static charges, one might hope to see the unstable string as an excited 
'resonant' string in the string mass spectrum. 

In addition to the complete string breaking described above, the gluons may propa- 
gate only some short distance along the string before returning and annihilating. These 
virtual processes will renormalise a-ji, and simple theoretical expectations for the string 
tension need to take this effect into account. 

Since we are considering larger SU(A^) groups (partly in order to calculate ak for 
larger k) one immediate question is how this screening will depend on A^. In particular 
we know that particle decay widths vanish in the large- A^ limit 0] and so it is natural 
to ask if screening will vanish in a similar way. The answer is yes and no. To appreciate 
this consider, say, the decay p — > 27r in large-A^ QCD. This is suppressed by a factor 
of 1/A^. However this suppression does not arise from the decay per se, but is a 
consequence of confinement constraining the pions to be colour singlets. If the theory 
were not confining, so that the 'vr'-mesons belonged to the adjoint representation of the 
colour group, then this decay of the p would be unsuppressed once we summed over all 
the coloured 27r final states. Thus the large- A^ suppression of particle decays can be 
thought of as a phase-space suppression due to confinement. In just the same way the 
process of gluon screening (and renormalisation) of strings will be unsuppressed at large 
A^. However the screening of a string in representation 7?. to a particular representation 
TZ' in the same AA-ality class may be suppressed. Whether it is or is not will depend on 
the number of states in 7^'. So, for example, adjoint string breaking, i.e. the adjoint 
sources being screened by gluons to singlets, will be suppressed as A^ — > oo. So will be 
the screening of fc = 1 strings down to the fundamental and in general the screening of 
/c-strings to the representation with k quarks. On the other hand the transformation 
of the mixed to totally antisymmetric k = 2 representations is not suppressed. (See 
Appendix [A.2| for details.) Of course these general counting arguments should be 
supplemented by any dynamical information we have. For example we expect ak —>■ ka 
as N —^ oo, from the suppression of fluctuations in that limit (and the dominance of a 
single Master field). This has implications for decays as well. 

3 Lattice preliminaries 

The way we perform our lattice calculations is entirely standard and follows the pattern 
described in 0. For completeness we shall provide a brief summary here. 

We shall work on a hypercubic lattice with periodic boundary conditions. The 
degrees of freedom are SU(A^) matrices, Ui, residing on the links, /, of the lattice. In 
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the partition function the fields are weighted with expl^} where S is the standard 
plaquette action 

S = -PY.{^-j^ReTrU,y (7) 

i.e. Up is the ordered product of the matrices on the boundary of the plaquette p. For 
smooth fields this action reduces to the usual continuum action with /3 = IN/ g"^ in 
D=3+l and (3 = 2N/ag'^ in D=2+l (where g^ has dimensions of mass and the theory 
is super- renormalisable). By varying the inverse lattice coupling (3 we vary the lattice 
spacing a. 

The Monte Carlo we use mixes standard heat-bath and over-relaxation steps in the 
ratio 1 : 4. These are implemented by updating SU(2) subgroups using the Cabibbo- 
Marinari prescription |ll|. We use 3 subgroups in the case of SU(3), 6 for SU(4), 10 for 
SU(5) and 15 for SU(6). To check that we have enough subgroups for efficient ergodicity 
we use the same algorithm to minimise the action. We find that with the above number 
of subgroups, the SU(iV) lattice action decreases more-or-less as effectively as it does 
in the SU(2) gauge theory. We calculate correlation functions every 5'th sweep. 

We calculate correlations of gauge-invariant operators 0(t), which depend on field 
variables within a given time-slice, t. The basic component of such an operator will 
typically be the (traced) ordered product of the Ui matrices around some closed contour 
c. A contractible contour, such as the plaquette itself, is used for glueball operators. 
If, on the other hand, we use a non-contractible closed contour, which winds once 
around the spatial hyper-torus, then the operator will project onto winding strings of 
fundamental fiux. In the confining phase the theory is invariant under a class of centre 
gauge transformations and this ensures that the overlap between contractible and non- 
contractible operators is exactly zero, i.e. the string cannot break. For our lattice 
action the correlation function of such an operator has good positivity properties, i.e. 
we can write 

C{t) = (0t(t)0(O)) = Y: I {n\<l>\n)\'exp{-Ej} (8) 

n 

where \n) are the energy eigenstates, with En the corresponding energies, and 1^2) is the 
vacuum state. If the operator has (0) = then the vacuum will not contribute to this 
sum and we can extract the mass of the lightest state with the quantum numbers of 0, 
from the large-t exponential decay of C{t). To make the mass calculation more efficient 
we use operators with definite momentum. (We will often use p = 0; however, as we 
will see, when better precision is required, it can be useful to extract extra information 
from the smallest non-zero momenta.) Note that on a lattice of spacing a we will have 
t = arit, where nt is an integer labelling the time-slices, so that what we actually obtain 
from eqndl) is aEn, the energy in lattice units. 

In practice a calculation using the simplest lattice string operator is inefficient be- 
cause the overlap onto the lightest string state is small and so one has to go to large 
values of t before the contribution of excited states has died away; and at large t the 



signal has disappeared into the statistical noise. There are standard methods |0 for 
curing this problem, using blocked (smeared) link operators and variational techniques. 
Here we use the simple blocking technique described in detail in 0]. We then have a 
set of trial operators corresponding to different blocking levels. From the space of oper- 
ators spanned by these we can determine the best operator using standard variational 
techniques |]1|. 

Having determined our 'best' operator, we then attempt to fit the corresponding 
correlation function, normalised so that C{t = 0) = 1, with a single exponential in 
t. (Actually a cosh to take into account the temporal periodicity.) We choose fitting 
intervals [^1,^2] where initially ti is chosen to be ti = and then is increased until an 
acceptable fit is achieved. The value of ^2 is chosen so that there are at least 3, and 
preferably 4, values of t being fitted. (Since our fitting function has two parameters.) 
Where ti = and the errors on C{t = a) are much smaller than the errors at t > 2a, 
this procedure provides no significant evidence for the validity of the exponential fit, 
and so we use the much larger error from C{t = 2a) rather than C{t = a). (This 
typically only arises on the coarsest lattices and/or for very massive states.) We ignore 
correlations between statistical errors at different t and attempt to compensate for 
this both by demanding a lower x^ for fhe best acceptable fit and by not extending 
unnecessarily the fitting range. (Although in practice the error on the best fit increases 
as we increase the fitting range, presumably because the correlation in t of the errors is 
modest and the decorrelation of the operator correlations is less efficient as t increases.) 
The relatively rough temporal discretisation of a few of our calculations, means that, 
at the margins, there are inevitable ambiguities in this procedure. These however 
decrease as a -^ 0. Once a fitting range is chosen, the error on the mass is obtained by 
a jack-knife procedure which deals correctly with any error correlations as long as the 
binned data are statistically independent. Typically we take 50 bins, each involving 
between 2000 and 40000 sweeps depending on the calculation. It is plausible that bins 
of this size are independent; however we have not stored our results in a sufficiently 
differential form that we can calculate the autocorrelation functions so as to test this 
in detail. A crude test is provided by recalculating the statistical errors using bins that 
are twice as large. We find the errors are essentially unchanged when we do so, which 
provides some evidence for the statistical independence of our original bins. 

In addition to the tension of the fundamental k = 1 string we also calculate tensions 
of fc = 2 and k = 3 strings. Denote by Pc the ordered product of the Ui around a 
non-contractible loop c that winds once around the spatial torus. So TrPc will project 
onto a winding loop of fundamental flux. The operators TrP^ and {TrPc}^ will project 
onto k = 2 loops, while the operators TrPj?, TrPc{TrPc}^ and {TrPc}^ will project 
onto k = 3 loops. These operators together with the same ones using blocked links, 
are summed so as to have p = and are then used as the basis of our variational 
calculation for the k = 2 and k = 3 strings respectively. 

We shall frequently perform fits to our extracted masses or string tension ratios. 



This may be to extract a continuum limit, or to calculate the Liischer string correction. 
The best fit is obtained by minimising the x^ and the value of the latter is used to 
determine whether the fit is acceptable. Equally conventional is our estimate of the 
error. Suppose we wish to calculate some quantity A through the fit. Let us suppose 
that the calculated 'data' values are Xj with cerresponding errors ai. An estimate of 
the error a a is provided by 

When the errors are small this formula is adequate, but may become unreliable for 
larger errors and poorer fits. It is however widely used and we therefore adopt it to 
facilitate comparison with other work. 

4 A universal string correction? 

Whether the long-distance dynamics of a confining flux tube is described by an effec- 
tive string theory and, if so, what is its universality class are fundamental theoretical 
questions. These are also important practical questions; particularly for an accurate 
determination of the string tension, since the answer will determine how large is the 
slowly falling 0(1//) correction to the mass of a long flux tube in eqns(3,|D. These are, 
however, difficult questions to answer numerically requiring, as they do, the accurate 
calculation of flux tube masses when these are very long and very massive. Thus, 
although this problem has been addressed many times in the past, the numerical evi- 
dence is, as yet, far from convincing. In this Section we shall describe some calculations 
which aim to improve significantly upon this unsatisfactory situation. 

Ideally we would like to perform calculations for the various SU(A^) gauge groups 
that are of interest to us in this paper. In practice our limited computational resources 
force us to focus upon the SU(2) group. Ideally, again, we would wish to perform 
calculations for several values of a but again this is not practical. Instead we shall 
perform calculations at a single value of a which is small enough, a^Ja ~ 0.16, that we 
can be confident that we are on the weak-coupling side of any roughening transition. 
We shall perform such calculations separately for 2+1 and 3+1 dimensions since both 
cases are of interest and they need not be the same. 

When is a string 'long'? Since we expect S,a = Wcr to provide the natural length 
scale for the physics of the confining flux tube, a string of length I = aL will be long if 
l/^cr = Lay/a S> 1. We can translate to more familiar physical units by recalling that in 
the real world l/\/cr ^ 0.45fm. Since quenched QCD provides a good approximation to 
QCD, we can, for qualitative purposes, use the same scale in the D=3+l SU(3) gauge 
theory. Since it appears that all D=3+l SU(A^) gauge theories are 'close' to each other 
P, it should not be too misleading to use the same scale in all cases. For purposes of 
orientation (and nothing else) we shall also use this scale in D = 2 + 1 SU(A^) gauge 
theories (where again all SU(A^) gauge theories are 'close' to each other |T|). 
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4.1 SU(2) in 2+1 dimensions 

We perform calculations on LF'Lt lattices at /3 = 9.0. (Recall that in D = 2 + 1 
SU(A^) gauge theories fj] the coupling g"^ has dimensions of mass, the theory is super- 
renormalisible, and j3 — »• 2N/ag'^ as a —>■ 0.) The flux tube winds around the spatial 
torus and so has length / = aL. We perform calculations for a large number of lattice 
sizes, ranging from L = 8 to L = 40. Recall that at this value of /3 one has a^/a ~ 0.162 
[|l| so that the length of our flux tube ranges from / = 8a ~ 1.3 x ^o- to / = 40a ~ 6.5 x ^o-. 
In our 'fermi' units the latter translates to / ~ 3fm. This should certainly be long 
enough to be governed by the long distance effective string dynamics if that indeed 
provides the correct description. 

In addition to the masses, as extracted from the p = operators, we also calculate 
the energies corresponding to the lowest five non-zero momenta transverse to the string: 
ap = 2nn/L for n = 1, .., 5. If we want to use these energies to provide extra informa- 
tion on the flux loop mass, care is needed because the continuum energy-momentum 
dispersion relation, E"^ = p"^ + m^, suffers lattice corrections. To determine these we 
have fitted the energies on our largest lattices with a more general dispersion relation 

{aEY = {amf + {apf + j{ap)^. (10) 

We find that for the largest lattices the size of the lattice correction 7 is consistent 
with zero within small errors. For example on the L = 40 lattice we obtain 7 = 0.06(8) 
for n < 5. Thus in these cases we can simply set 7 = and use the continuum 
dispersion relation. For smaller lattices the gap between momenta becomes larger and, 
not surprisingly, the number of momenta that can be well fitted with E"^ = p'^ + m^ 
decreases. Since these larger values of E have larger statistical errors there is not much 
to be gained by attempting to include them and so we simply exclude them from the 
fits. For the same reason we do not bother with p 7^ for lattices with L < 20. 

We note that the lattice correction in eqn(|TDp is 7(ap)^ = '^{2im/LY and this is of 
the same order as the higher order non-universal string corrections that we have ignored 
in eqn(|D. We also note that the tree-level lattice dispersion relation provides very bad 
fits to our calculated energies - it clearly contains lattice corrections that are far too 
large. Finally we remark that if we generalise eqn ([IO|) so as to allow a renormalisation 
of the O(p^) term, we find that the fitted coefficient is unity within very small errors. 

We list the masses that we obtain from our p = correlators in Table |l|. To obtain 
the values in the first column we have used fits to the correlation functions, C(t), down 
to the lowest plausible values of t, so as to minimise the errors (which grow with t). 
In some cases there are indications from the effective masses at larger t that this mass 
estimate might be optimistic. So we have listed in the second column mass estimates 
that we regard as very safe, but which might, as a result of being overcautious, overstate 
the errors and hence weaken the statistical significance of our final fits. Table |l| also 
contains, in the third column, the mass estimates obtained by using both p = and 
p 7^ energies in the way described above. Statistically these are the most accurate 
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results, although they run the risk of possessing a small systematic bias from lattice 
corrections to the continuum energy-momentum dispersion relation. However any such 
bias will be smallest on the largest lattices and it is only on these that the p ^ 
energies make a significant difference. 

An immediately striking feature of the listed masses is that they rise more-or-less 
linearly with length, all the way out to the longest loops. This demonstrates directly 
that in D=2+l SU(2) gauge theories we have linear confinement out to at least / ~ 3fm. 
However this is no more than expected and so we shall not dwell upon it any further 
here. 

We turn now to the real question of interest here: how accurately can we test the 
0(1/1) string correction term in eqn(|^)? As a first step we calculate the effective 
value of the coefficient c^ that one obtains from pairs of flux loops of length / and /' 
respectively: 

^ r m{l) rail') t 

c:^^(MO = -Tyf^ X ^7; - 7-^ (11) 

In Table @ we list the values of dg^^il^l') that one obtains for neighbouring values of 
I and /' using the masses listed in Table |I]. Any given value gives us no information 
on the validity of the 0{l/l) string correction. However if we find that c1^^{l^ V) has a 
finite non-zero limit as Z, Z' — »• cxo then we will have shown that the leading correction 
is indeed of this form and the limiting value will provide us with an estimate of the 
coefficient c^. 

Our most accurate values of d^J^ are those that incorporate p 7^ energies, and we 
display these in Fig.|I]. We see from the plot that for small loop lengths the value of 
dl/f is small and increases as the loop length grows. However the behaviour is not 
monotonic: at intermediate / the value of c^J^ increases through the bosonic string 
value and perhaps peaks close to the value for a Neveu-Schwartz string. This occurs 
at a loop size / ~ Ifm which is typically the longest loop for which older calculations 
had usefully accurate results. Thus, by focusing on slightly different intervals close to 
I ~ Ifm it is possible to either confirm or contradict the bosonic string value; but, in 
either case, incorrectly. To obtain real evidence one must go to longer strings and if 
one does so, as in Fig.|I], one finds that the value of c^J^ decreases again to something 
that appears consistent with the bosonic string value. 

To obtain our estimate of Cs we use eqn(^ to fit all the loop masses that are longer 
than some reference value /q- If ^0 is small there is no acceptable fit. As we increase /q 
eventually the fit becomes acceptable. We can then increase /q to check the stability of 
the best fit. In Table |^ we list the results of this procedure for each of our three sets of 
loop masses. The most accurate values are obtained from the last set, and are plotted 
in Fig.^]. From the second column of Table |^ we extract a 'safe' estimate for the string 
correction: 

c, = 1.066 ± 0.036. (12) 

This is close to the bosonic string value and far from that of other simple string theories. 
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The above analysis assumes that the leading correction to the linear dependence 
of the mass is oc l/l. The fact that d'J^ becomes independent, within errors, of I 
once / > 16a, tells us that our results are certainly consistent with such a string-like 
correction. It is interesting to ask how well our calculations exclude other choices. We 
have therefore performed fits to 

m(0 = a/ + |. (13) 

We find that the there are no acceptable fits if any values with / < 14a are included. 
The fit to Z > 14a has a mediocre but acceptable x^) and we find that the possible 
range of powers is p = 1.4 ± 0.5. For I > 16a the best fit is very good and one finds 
p = 0.9 lb 0.5. Fits to / > 20a are equally good but no longer provide a useful constraint 
on p. In short, our results are consistent with the 0(1//) string-like correction term 
and, in any case, the power of l/l is constrained to be within the range p = 0.9 ± 0.5. 
So if we constrain the power p to be an integer, we find that it indeed has to be p = 1. 

4.2 SU(2) in 3+1 dimensions 

We perform calculations on L^Lt lattices at /3 = 2.55. At this value of (3 one has 
a^/a ~ 0.159 0, so that the size of a is very similar to that in our D = 2 + 1 
calculations. The calculations in D=3+l are, of course, slower and so our range of 
lattice sizes and our statistics is somewhat less. One may hope that this will be partly 
compensated for by the fact that the expected string correction, {D — 2)ttCs/QL, will 
be twice as large (for a given universality class). We perform calculations on lattices 
ranging from L = 8 to L = 32. Thus our longest flux loop is / = 32a ~ 2.3fm which, 
if our experience in D = 2 + 1 is relevant, should be long enough to see the leading 
correction. 

Our calculation and analysis is precisely as in D=2+l, except that the values of the 
momenta transverse to the flux loop that we use are p"^ = 0,1,2,4. The mass estimates 
are listed in Table ^. As in the D=2-|-l case we list two sets of masses extracted from 
the p = correlators. In general our mass estimates are chosen to be those with the 
smallest errors while still giving plausible fits. In some cases the plausibility is less 
than convincing and we then also select a 'safer' mass estimate, which will have larger 
errors. The former numbers provide the first column of masses in Table § while the 
latter provide the second column. The two columns only differ in some cases. The mass 
estimates obtained using p ^ as well as p = correlators are also divided into two 
sets. (A division that did not appear useful in D=2+l.) The first set (third column of 
masses) differs from the second principally in that on the L = 14, 16 lattices we chose 
less plausible p = masses in order to be consistent with the p = 1 values with which 
they were then averaged. In the last 'safe' column we dealt with this discrepancy by 
not using the p ^ values (which, in any case, become much less useful on the smaller 
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lattices). Thus this range of mass estimates gives some indication of any systematic 
error that arises from our procedure for extracting masses. 

We first note that the loop mass increases approximately linearly with the loop length 
confirming, as expected, that the theory is linearly confining (up to ~ 2.3fm). 

In Table ^ we list the values of the effective string correction coefficient, c^-^-^, defined 
in eqn(^). As in D=2+l our most accurate values of c\^^ are those that incorporate 
p 7^ energies, and we display one set of these in Fig.|^. We see in the plot a behaviour 
similar to what we observed in D=2+l: the value of d^J^ increases as the loop length 
increases, attains a maximum value at / ~ Ifm that is significantly larger than the 
bosonic string value, and then decreases to a value consistent with the value for a 
bosonic string. 

Just as in D=2+l we estimate Cg by using eqn(^ to fit all the loop masses that are 
longer than some reference value Zq- If ^o is small there is no acceptable fit. As we 
increase Iq we find that the fit eventually becomes acceptable. We then increase Iq to 
check the stability of the best fit. In Table ^ we list the results of this procedure for 
each of our last three sets of loop masses. (For the first set there are no acceptable 
fits, perhaps indicating that some of the mass choices were indeed too optimistic.) The 
most accurate values are obtained from the last two sets, and we use these to obtain 
our best estimate for the string correction: 

c, = 0.94 ±0.04. (14) 

This is close to the bosonic string value and far from that of other simple string theories. 

5 /c-strings 

In the previous Section we accumulated some evidence that SU(2) flux tubes in the 
fundamental representation are described by an effective bosonic string theory at large 
distances. In this Section we consider flux tubes in higher representations: the fc-strings 
described in the Introduction. We will not be able to perform comparable checks on the 
stringy nature of these flux tubes although we will perform some crude finite volume 
analyses that are primarily designed to confirm the presence of linear confinement in 
SU(iV > 4) gauge theories. In our analysis of the string tension ratios, cr^/cr, we shall 
make the plausible assumption that the leading correction is that of a simple bosonic 
string. However for completeness we shall also point out how the results are affected 
if this should not be the case. 

In order to provide some theoretical context within which to view our numerical 
results, we shall first briefly summarise some of the existing ideas about how such ratios 
might behave. This is not intended to be an exhaustive review, and our references 
are merely designed to provide an entry into the literature rather than aiming at 
completeness. 



13 



We then describe our calculations oi k = 2 and (fundamental) k = 1 strings in 
D=3+l SU(4) and SU(5) gauge theories. We follow this with a description of our 
D=2+l calculations which are for SU(4) and SU(6). In this last case we also have 
non-trivial k = 3 strings that we are able to analyse. 

We shall find that in D=3+l the string tension ratios lie between the predictions 
of MQCD and Casimir scaling, straddling both within two standard deviations; with 
SU(4) slightly favouring Casimir scaling and SU(5) leaning towards MQCD. In D=2+l 
our results are again close to both Casimir scaling and to MQCD, but now they are 
much more accurate and so we can begin to see significant deviations. Although we 
see deviations from both sets of predictions, those from MQCD are much larger than 
those from Casimir scaling. We point out that near-Casimir scaling occurs naturally 
if the confining flux tube has a cross-section that is nearly independent of the flux 
carried. We perform explicit calculations in D=2+l that suggest that this is in fact 
so. These calculations give us, as a side-product, the value of the /c-string tensions at 
high temperature, and we find near-Casimir scaling there as well. Motivated by this 
we perform a high T calculation in the D=3+l SU(4) gauge theory where we again 
find near-Casimir scaling. 

It is interesting to ask if all this also occurs for the unstable strings. We shall show 
that our D=2+l calculations provide some evidence that points to this . 

5.1 Some expectations for k-strings 

The interest in strings that emanate from sources in higher representations goes back a 
long way. The early discussions were framed in terms of unstable strings in SU(2) and 
SU(3) gauge theories as were the lattice calculations. (See, for example, 0.) Despite 
the uncertainties of this kind of calculation, these early results were already seen as 
being able to discriminate against particular theoretical ideas; in particular |T3| against 



the bag model [13, 14|. There have been recent much more accurate SU(3) calculations 
|T5| , P^ that support this earlier work, and this has sparked some interest in the possible 
dynamics Jl^ [18|, |1^. The recent interest [|20|, ||, including our own, in stable fc-strings 



in SU(iV > 4) has been largely due to conjectures arising in M(-theory)QCD [0]. Here 
we briefly allude to some of these theoretical ideas, with a particular focus on MQCD 
and 'Casimir scaling' since their predictions turn out to be closest to the results we 
obtain p, ^ for the tensions of fc-strings. 

5.1.1 unbound strings 

The simplest possibility is that 

ak = kak=i ; k = mm{k,N — k}. (15) 

That is to say the total flux is carried by k (or N — k ii that is smaller) independent 
fundamental flux tubes. This would occur if the interaction between fundamental flux 
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tubes was so weak that there were no bound or resonant multi-string states. One may 
regard this as the trivial scenario with which the actual calculated values of au can be 
contrasted. 

5.1.2 Casimir scaling 

The idea that the confining fiux tube between sources will be proportional to the 
quadratic Casimir, C-ji, of the representation of those sources 

a-n oc Cn (16) 

is an old idea. An early motivation arose from a model of 'random fiuxes' for 
the vacuum and the observation that in certain solid state systems this leads to a 
dimensional reduction D -^ D — 2. Thus D = 4 theories would reduce to D = 2 
gauge theories in which the Coulomb linear potential is indeed proportional to the 
quadratic Casimir. The numerical calculations supporting this were in D=4 SU(2) [Q 
and involved potentials at relatively short distances. The observation soon after that 
one seemed to see a similar Casimir scaling in D=3 theories 0] forced a generalisation 
of the dimensional reduction idea 0. The most accurate early calculations involved the 
adjoint string tension. Recently however there have been accurate calculations IT^, ITB 



for a variety of representations in SU(3) and this has sparked renewed interest in this 



idea [O, M 19 



There are obvious ambiguities in calculating the string tension of unstable strings 
from the intermediate distance behaviour of the static potential. At short distances 
we know that we have a Coulomb potential which, of course, displays Casimir scaling. 
As the potential interpolates between this and the long-distance behaviour one expects 
some continuity. If, as is usually done, one fits the potential V{r) by a simple sum of 
a Coulomb term and a linear piece, V{r) = Vq + c/r + ar, and then performs the fit 
in a limited range of r immediately beyond the Coulomb region, then it might be that 
simple continuity artificially forces approximate Casimir scaling on the fitted linear 
term. While this is no more than a possibility, it does underscore the utility of using 
stable strings, as we shall do, where one can go to larger distances, and doing the 
calculation in a way, as we shall also do, that does not involve explicit sources and 
associated Coulomb terms. 

Since the sources may be screened by gluons, which are in the adjoint representation 
and do not feel centre gauge transformations, it is appropriate, as we remarked earlier, 
to categorise the representations of SU(A^) sources by how they transform under the 
centre of the group. If the source acquires a factor z'^ with k G Zn, then we shall 
generically refer to the corresponding fiux tube as a A;-string. Any fc-string can be 
transformed into any other fc-string by appropriate gluon screening. Thus the stable 
/c-string will be the one with the smallest string tension. Any other fc-string will, at 
sufficiently large distances, find it energetically favourable to be transformed into the 
lightest stable string through gluon screening. If we have Casimir scaling then the 
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lowest string tension corresponds to the representation with the smallest quadratic 
Casimir, and this is the totally antisymmetric representation. (See the Appendix [A.3| .) 
The ratio of these quadratic Casimirs gives us the Casimir scaling prediction for stable 

/c-strings 

a, _ kiN - k)_ ^^^^ 



a N 



5.1.3 MQCD 



A number of calculations in brane (M-)theory of QCD-like theories (see and refer- 
ences therein), which are generically referred to as MQCD, find that that the string 
tension of fc-strings satisfies 



o-fc _ sm ^ 



kn 

(li 



a sm^ 

This led to the conjecture P| that this might be a universal result and that QCD (and 
SU(iV) gauge theories) fall into this universality class. 

This prediction has reasonable properties: it has the required k ^^ N — k symmetry 
and takes sensible values for A^ = 2, 3. However the MQCD derivation neglects po- 
tentially important quantum fluctuations which might 0] renormalise the simple and 
elegant formula in eqn(|T8|). 

The MQCD calculations are, strictly speaking, for SU(A^) gauge theories in 3+1 di- 
mensions. It is not clear how much evidence there would be for a corresponding MQCD 
conjecture in D=2+l, although a naive reading suggests that the brane constructions 
in [0 would lead to the same conclusion for the ak ratios. In any case the trigonometric 
formula in eqn(|l^) has the correct qualitative properties and so we shall compare our 
results to it not only in D=3+l but also in D=2-t-l and, indeed, at finite temperature. 

5.1.4 bag model 

In the bag model (see e.g. |]T1|, |T^) the flux between distant sources is confined to a 
cylindrical bag of cross-section A. The flux is homogeneous, 

EaA = gTa, (19) 

and the vacuum energy difference between the inside and outside of the bag is given 
by the bag constant B. Thus the energy per unit length is n, [131 



- = 2Tias-^ + AB (20) 

where C-ji is the quadratic Casimir of the source and as is the strong coupling constant. 
One now fixes the area A by minimising the energy. This gives the string tension to be 

cTn « {CnV^ (21) 
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which differs markedly from Casimir scahng. The fact that the early numerical calcu- 
lations gave an adjoint string tension that satisfied eqn(|16D rather than eqn(|2T|) was 



picked up [jT3| as providing critical evidence against conventional bag dynamics [|1J, , 
in that it suggested a flux tube cross-section that was independent of the size of the 
flux. 

5.1.5 strong coupling 

In the strong coupling limit, /3 — > 0, of our action, a Wilson loop involving k strings 
will need to be tiled with plaquettes at least k (or N — k) times. The leading term in 



this limit will therefore reproduce eqn(|l^): a^ = kak=i, k = nim{k,N — k}. However 
the non-leading terms will introduce interactions between these tiled surfaces, and this 
simple ratio will change as we move away from the strong coupling limit. 

Strong coupling predictions are, of course, not universal; however this one is more 
universal than most. If we generalise the action to contain any combination of closed 
loops, so long as these are linear in the SU(iV) link matrices we will still obtain eqn(|T5|). 
However if we include loops or products of loops that are not linear in the links then 
we can obtain other results. One can think of the action as having loops in different 
representations, and the value of Ck/cr will depend only on what these representations 
are and what are their relative weights. By choosing an action in an appropriate 
'universality' class, one can essentially obtain for cTfc/cr any value one wants. 



Hamiltonian strong coupling (see e.g. |21|) is more interesting. The leading term, 
as (7^ — s> oo, is simply the quadratic Casimir for each spatial link. Gauss's law means 
that our two fc-sources are joined by excited links, and that the lightest k string will 
satisfy Casimir scaling as in eqn(^). Of course the magnetic perturbation will spoil 
this result as we move away from /3 = 0. 

5.2 /c-strings in D=3+l 

We will now calculate the ratio of the k = 2 and fundamental string tensions, ak=2/o', 
in both SU(4) and SU(5) gauge theories. There are no other stable /c-strings for these 
values of N but having results for two values of A^ will already provide significant 
constraints. We are, of course, interested in the continuum limit, so we calculate this 
ratio for several lattice spacings and then extrapolate to the continuum limit using the 
fact that for the plaquette action the leading lattice correction to dimensionless mass 
ratios is 0{a'^): 

—— = — — + ca a. 22 

cr(a) a(0) 

We calculate the string tension from the mass of a flux loop that winds around the 
spatial torus. We assume that the leading correction to the linear dependence of the 
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mass is that appropriate to a simple bosonic string: 

mk{l) = CTkl y. (23) 

This assumption has some support from the calculations of the previous Section, but 
it is not guaranteed that what holds for SU(2) holds also for larger A^. So we shall 
occasionally pause to state how sensitive are our results to this assumption. 

We begin by listing in Table |^ the (fundamental) string tensions corresponding 
to the various (3 values at which we perform our calculations. This sets the scale of a 
in physical units. In Table |^ we list our lattices and calculated values of the k = 1 and 
k = 2 flux loop masses for the case of SU(4); and in Table ^ for SU(5). (Note that in 
these calculations we only use p = correlators.) 

In order to extract a string tension from the flux loop mass we must ensure that our 
loop length is long enough for the corrections to eqn([23|) to be negligible within our 
statistical errors. In Section K2 we have seen that in the case of SU(2) this appears 



to be the case for strings longer than ly/a = La^/a ~ 3 (see Tables ^^. Here we 
perform an additional finite size study, this time in SU(4), which, while less accurate, 
will probe the behaviour of fc = 2 as well as k = 1 strings. 

Our finite size study is at /? = 10.7 and involves loops ranging from L = 6 to L = 16 
with masses as listed in Table ^. The longest length translates into / ~ 4.9/^/o' ^ 2.2fm. 
We observe that both the k = 2 and k = 1 masses grow approximately linearly with /, 
demonstrating that the SU(4) theory linearly confines both k = 1 and k = 2 charges 
(at least over this distance range). Using eqn(|23|) we extract the ratio (7^=2/0" which we 
plot in Fig.|^. We see that within errors the ratio becomes independent of the flux loop 
length for I > 10a c^ 3/-\/o". For comparison we also show what happens if we do not 
include any string correction at all, i.e. 17^=2/0" = 771^=2 (0/'"^fc=i(0- We see that while 
the ratio changes by a few percent, it becomes independent of /, within errors, at the 
same length, / = 10a. By the same token it is clear that the results of this calculation 
are not accurate enough to distinguish between different possible string corrections. 

Our finite volume study has taught us that higher order corrections in 1/Z to the 
string tension ratio will be negligible (within our typical errors) if we make sure that 
our loop length satisfies /a/c" > 3. Comparing the values of a^/a in Table |^ with the 
corresponding lattice sizes listed in Table ^ and Table |^ we see that our loop lengths 
have been chosen to fulfill this bound; more generously at smaller a where the errors 
are smaller. 

Assuming eqn(p3D, we extract our string tension ratios and plot them against a'^a in 
Fig. 1^. (At f3 = 10.7 we use only the L = 10, 12 lattices since the larger volumes have 
errors that are too large to be useful.) On such a plot the continuum extrapolation, 
eqn(p2D, is a simple straight line. We show the best such fits in Fig. |^. We find that if 
we use all the points we get an excellent x^ fo^' SU(4) and an acceptable one for SU(5). 
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From these fits we obtain the following continuum values: 

a^_r 1.357 ±0.029 SU(4) 
a™ a \ 1.583 ±0.074 SU(5) ' ^ > 

One might worry that these fits could be biased by including the coarsest a value 
(where we know the lattice corrections to the scalar glueball mass to be large). If 
we exclude this coarsest a point our best values in eqn(P3) are changed to: 



a^=2 _ J 1.377 ±0.035 /?> 10.70 SU(4) 
i™ a ~\ 1.76 ±0.14 /5> 16.975 SU(5) ' ^^ 

This gives us some idea of the direction of any such bias. 

It is interesting to compare our results with the expectations of MQCD 



and Casimir scaling 



(Tk=2 MQCD Sin-^ 


1.41... 


SU(4) 


a sin;| 


1.61... 


SU(5) 


ak=2 cs k{N - k) 


j 1.33 


SU(4) 


a N-l 


1 1.50 


SU(5) 



(26) 



(27) 



We see that our results in eqn(|21) and eqn(P^) are consistent, at the 2a level, with 



both these expectations, within quite small errors; with perhaps a slight bias towards 
favouring MQCD. It is because the two sets of predictions are numerically very similar 
that we cannot, at present, choose between them. On the other hand we clearly exclude 
the unbound string value of 2: i.e if we do wish to think of the k = 2 string as being 
composed of two k = 1 strings then it must be a tightly bound state of such strings. 
We also clearly exclude the bag model prediction: 



(^k=2 Bag j kjN - k) _ f 1.15... SU(4) 

a ~ V iV-1 ~ \ 1.22... SU(5) ' ^^^^ 

In order to distinguish clearly between MQCD and Casimir scaling we need to reduce 
our statistical errors by about a factor of two; a feasible goal but one for the future. 

Thus our conclusions are essentially unchanged from those of our earlier paper [Q 
although our SU(4) calculation now has smaller statistical errors, and our SU(5) cal- 
culation is now free of the potentially large systematic errors that concerned us earlier. 

A final remark. Our above analysis assumed that the flux tubes behave like simple 
bosonic strings. What difference does it make if we do not make this assumption? 



Suppose first that we use the result we obtained in Section |4]^ for the coefficient of 
the string correction: Cg = (1.25 ± 0.25)7r/3 (where we take very generous errors). 
Repeating our analysis with such a string correction we find that our results in eqn(p^) 
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are lowered by about 10% of the statistical error; that is to say, insignificantly. Even if 
we were to ignore what we knew and simply assumed some range like c^ = (1 ± l)7r/3 
we would find that the maximum shift would be less than our quoted statistical error. 
(For Cg = the ratios are close to MQCD while for Cg = 2tt/3 they drop very close to 
Casimir scaling.) 

5.3 A;-strings in D=2+l 

Our calculations in D=2+l follow the same pattern as in D=3+l except that our 
calculations are in SU(4) and SU(6). The main reason for SU(6) rather than SU(5) is 
that with the former one also has stable k=3 strings that one can study. On the other 
hand the calculations take longer which is why we contented ourselves with SU(5) in 



four dimensions. Our calculations are summarised in Tables ^ and |12] and the scale 
of a, in units of the fundamental string tension, is given in Table |10|. 

We begin with a finite volume study in SU(4) at /5 = 28.0 that parallels our D=3+l 
study. The loop lengths range from L = 4 to L = 16, with the largest loop correspond- 
ing to I = 16a ~ ^1 ^(y — 1.6fm. We extract ak=2/o' using eqn (|2^) and plot the result 
in Fig.^. We see that the ratio of string tensions is independent of the loop length 
(within errors) once / > 10a ~ l.'bj \ja ~ l.lfm. This is a somewhat shorter length 
than the one we found in D=3-|-l. We shall later see that the flux tube is thinner (in 
units of 0") in D=2+l than in D=3+l and this is presumably why the corrections are 
smaller. We also show in Fig.^ the string tension ratios one obtains if one assumes no 
correction. This also plateaus for / > 10a. Moreover we see that the value of the ratio 
differs by only about 1%. (Note the string correction is oc (-D — 2) and so is larger in 
D=3+l than in D=2+l.) 

We observe that the loop lengths we shall use, as listed in Tables |TT] and [T^, satisfy 
the above bound, / > 2.5/1/0", by a good margin. So assuming eqn(p3D, we plot our 
string tension ratios in Fig. |^ against a^a. We also show the straight line continuum 
extrapolations, using eqn(p2D. We find that we get an acceptable x^ using all the 
points. We thus obtain the following continuum values: 

o^ ^ r 1.3548 ± 0.0064 SU(4) 
a"o a \ 1.6160 ±0.0086 SU(6) ^ ^ 

and 

lim^^ = 1.808 ±0.025 SU(6).. (30) 

We note that the errors here are much smaller than in eqn(^) and the MQCD expec- 
tation is excluded. We also see, in our k = 2 SU(4) ratio, a deviation from Casimir 
scaling at the level of about 3.5 standard deviations. In SU(6) the expectations are 

ak=^sujp) f 1.73... MQCD 
a \ 1.60 CS ^ ^ 
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and 

ak=3 su{6) j 2.0 MQCD 



a I i.o CS • ^^^^ 

We see that in this case our resuhs are far from the MQCD values and agree well with 
Casimir scaling. 

Given our very small statistical errors one might worry that our assumption of a 
bosonic string correction might introduce a relatively large systematic error. In fact 



this is not so. If we take our SU(2) string analysis in Section ^?1| to be telling us that 
the coefficient of the l/l term is Cg = (1.1 ± 0.1)7r/6, then we find a negligible shift in 
our above predictions. Even if we were to assume c^ = this would shift our values 
upwards by less than 2%. While this would increase the deviation from Casimir scaling, 
it would be far from enough to bridge the gap to the MQCD prediction. However Cg = 
is a contrived and extreme example, and it appears to us that any reasonable estimate 
of the systematic error arising from the uncertainty in the string correction shows it to 
be negligible. 

5.4 Width of /c-strings 

We have seen that in both D=2+l and D=3+l the ratio cr^/cr is close to the Casimir 
scaling value. In D=3+l it is also consistent with MQCD, but the MQCD and Casimir 
scaling predictions are in fact quite close. If one imagines modelling the fiux tube then 
this is a somewhat conter- intuitive result. It would be natural to think of the fiux as 
homogeneous, as in eqn(|T^, and that the vacua inside and outside the fiux tube differ 
by some energy density 6E^. These are of course the ideas embodied in the Bag model. 
One would then expect that as the representation of the fiux increases, so that the 
chromoelectric energy density increases, the area will increase so as to minimise the 
total energy increase. This is just the variational calculation of the naive bag model 
which leads to a ratio cfk/cr that grows as the square root of the quadratic Casimir. 
This is definitely excluded by our calculated values. One might imagine extending this 
simple-minded model by providing the fiux tube with a surface tension. However this 
would have no effect in D=2+l where the 'surface' is independent of the flux tube 
width. 

If the flux is homogeneous and if the flux tube width is independent of the total flux 
carried, then one naturally obtains Casimir scaling. If one considers a superconductor, 
in a phase that exhibits the Meissner effect and supports (magnetic) flux carrying flux 
tubes, there is a range of parameter values, called the deep-London limit, where the flux 
tube cross-section is indeed independent of the flux (see e.g. [^). This corresponds to 
a penetration depth, related to the photon mass, being much larger than the inverse 
scalar Higgs mass. The deviation from Casimir scaling would be related to the ratio 
of these masses. 

It is interesting to test these ideas. Here we shall attempt to calculate the widths 
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of the flux tubes corresponding to different /c-strings and see how close the width is to 
being independent of k. We shall do so in D=2+l because it is faster; but the same 
technique can be used in D=3+l. We shall perform calculations for k = 1,2 flux tubes 
in SU(4) and for A; = 1, 2, 3 flux tubes in SU(6). 

We use a technique that was employed in pSf to calculate the width of SU(2) flux 
tubes in D=2+l. Consider a lattice of size L x L± x Lt. Lt refers to the Euclidean 
time in which we calculate correlations. The flux loop is of length L, and L± is the 
spatial size transverse to the flux tube. By reducing L± we can squeeze the flux tube. 
If the flux tube oscillates with simple harmonic modes then it will not be affected by 
reducing the finite (periodic) transverse width until it reaches the 'intrinsic' width of 



the flux tube [^]. Once L± is smaller than this width, which we shall call l^ = aL 



Wl 



we expect the mass of the flux loop to increase as ||23| 



am{L] L±) = ain[L; oo) x - — . (33) 

The onset of the increase is at L± = L^ and provides us with an estimate of L^. Our 
main interest is to see if Ly^ varies with k or not. 

Of course the above arguments are very simple. There are also changes in the vacuum 
as L± becomes small which we have neglected. In particular there is a phase transition 
at a critical value of L± (see Section |5.5| ) which is characterised by the development 
of a non-zero vacuum expectation value for the Polyakov loops that wind around the 
short L± torus. However the string tensions we calculate behave smoothly through this 
transition, suggesting that our simple analysis should not be invalidated. In any case, 
our main conclusion, that as k grows the flux tube width does not grow ever larger, 
will remain valid since such a growth would mean that higher k flux tubes would begin 
to be squeezed when L±_ was greater than its critical value, and this would certainly 
be visible. 



We show in Table ^ how amk{L;L±) varies with L± for k = 1 and k = 2 flux 
tubes; all at /5 = 28 in the SU(4) gauge theory in D=2+l. We do so for L_l = 2, ..., 20 
and for two values of the flux tube length, L. The minimum transverse size is l± = 
aL± = 2a ~ 0.5/^cr which we expect to be smaller than l^. The loop lengths are 
I = aL = 8a, 12a ~ 0.9, 1.35fm which should be long enough to allow well-formed flux 
tubes. The fact that we have calculated masses for two values of L at each L± allows 
us to check whether the mass is growing nearly linearly with L and that we do indeed 



have a flux tube. We see from the masses listed in Table ^ that this is so for all values 
of L^. 

We plot the L = 8 flux loop masses in Fig.^. We see that both mk=i and mk=2 
start increasing at a very similar value of L±. Moreover the masses for the smallest 
two values of L± are consistent with eqn(p3D. If we take these to fix the value of L 
in eqn (|33|) , we find that Ly^ ~ 5.2 for both k = 1 and k = 2 flux tubes. That is to say, 
the flux tube width is indeed independent of the flux, and has a value l^ ~ 1.3/^cr. 
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In addition to these gross features we also see in Table ^ and Fig.^ that there is 
a decrease in the mass at values of L± that are just above the values where the mass 
starts to increase. Moreover this decrease is more pronounced for the k = 2 string 
than for the k = 1 string. Such an effect indicates that there are some difference 
between the sizes of the two flux tubes - if only in their tails - and an analysis of this 
might provide information on the dynamics, e.g. the parameters of the effective dual 
superconductor referred to earher. However anything quantitative needs calculations 
with more resolution, i.e at smaller values of a. 

We display in Fig.^] how (7^=2/0" varies with L±. We see that the ratio is close to 
Casimir scaling not only at large L± (something we have seen already) but also at 
small L±. There is only a small range of L±, coinciding with the dip in mk, where 
the ratio drops below this. Of course, at very small values of L± our D=2+l system 
is effectively reduced to a D=l+1 gauge-scalar theory, and we recall that the linear 
confinement of pure gauge theories in D=l+1 arises through the Coulomb interaction 
which automatically satisfies Casimir scaling. 

Our SU(6) analysis closely parallels our SU(4) analysis except for the fact that we 
now have additional k = 3 strings. Our calculations are at /3 = 60 which, as we see 



from Table |T0|, has a similar a to that at /? = 28 in SU(4). 

Our masses are listed in Table |l^. They are for flux loops of length I = 10a and 
/ = 12a and we see evidence for a linearly growing mass for all k and for all L_j_. We 
plot the L = 10 flux loop masses in Fig.|lO| and we see, once again, an increasing loop 
mass at small L± that is consistent with eqn(^). Indeed we find a common flux tube 
width, Lw ~ 4.5, for all three values of k. So, just as in SU(4), the flux tube width is 
independent of the flux, and has a value l^ — 4.5a ~ \.1j ^o that is also very similar. 
Again, just as in SU(4), the loop mass decreases just before it begins to increase. 

We plot oi^ja in Fig.|TT[ We again see consistency with Casimir scaling at small as 
well as at large L^, except possibly in the region of the dip. 

This analysis thus appears to confirm that the confining flux tube has a cross-section 
that is largely independent of the flux carried. The minor differences between fc-strings 
might, however, be useful in telling us about the details of the confinement mechanism. 

5.5 High T spatial string tensions 

In our above calculations we have calculated the mass of a long flux loop in a spatial 
volume with a limited transverse spatial dimension, L^. Let us relabel the axes of our 
L X Lx X Lt lattice so that the short spatial torus becomes our time torus and our time 
torus becomes a long spatial torus. We are now on a L^. x Ly x L^ lattice with Lt = L±_ 
and Lx,Ly ^ Lj. This corresponds to a system at temperature aT = 1/L^. and, as 
we decrease L±, we will pass through the deconfining phase transition at T = Tc. In 
this rotated co-ordinate system the flux loop 'mass' that we have calculated is obtained 
from the spatial correlation of spatial loops; it is a screening mass, from which we can 
calculate what is usually referred to as the 'spatial string tension' (usually obtained 
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from spatial Wilson loops). Thus our finite width studies have in fact provided us with 
a calculation of the k = 1, k = 2 and k = 3 spatial string tensions as a function of T 



in SU(4) and SU(6) D=2+l gauge theories. (All this parallels previous studies |2^, ^ 
of SU(2) in D=2+l.) 

Simple arguments (see for example |2^) tell us to expect a oc g^T for T ^ Tc. In 
our case, where aT = 1/Lt = l/L±, this translates to a^a oc l/L± as L±_ — > 0. This 
is precisely what we have already inferred from the SU(4) and SU(6) calculations in 



Tables ^ and |14[ Indeed we see that the linear increase with T sets in very close to 
T = Tc, and does so simultaneously for all fc-strings - presumably due to the squeezing 
of a fiux tube whose width is ~ l/T^. (Although we do not have precise calculations 
of the deconfining temperatures in SU(4) and SU(6) gauge theories, we expect from 
extrapolations of previous SU(2) and SU(3) calculations that Tc ^ O.QS-y/cr in D=2+l 
and this tells us that the critical value of L± is ~ 4 at these values of (3.) Moreover, 



as we see from Figs. § and |1T], the string tension ratio is close to the value expected 
from Casimir scaling. This is not a great surprise: the high-T dimensional reduction 
of the D=2+l theory takes us to a D=H-1 theory, and in a D=l+1 gauge theory 
even the Coulomb potential is linearly confining; and the latter will automatically 
satisfy Casimir scaling. This is of course simplistic; the dimensional reduction leads 
to (adjoint) scalars as well as gauge fields, in the reduced theory and, in any case, the 
linear potential may have other sources than just the Coulomb potential. 

We now turn to the more interesting case of D=3+l. We shall not attempt a 
systematic study but, just as in D=2+l, we shall work at one single value of /3 and 
will vary T in the rather coarse steps allowed by varying Lt. Our calculation is in 
SU(4) at /5 = 10.7. Although we do not have precise information on the deconfining 
temperature, an extrapolation of previous SU(2) and SU(3) calculations |^ using a 
simple 0{1/N^) correction, suggests that Tc ^ 0.62^0" for SU(4). Since ay/a ~ 0.306 
at /? = 10.7 (see Table 0) the critical value of Lj is ~ 5. Thus a spatial volume of 10"^ 
should be adequately large for our exploratory calculation. Accordingly we work on the 
lO^Lt lattices listed in Table ^. In the Table we also include our earlier calculations 
on a 10^ lattice, to provide the 'T = 0' reference value. 

We first wish to establish the rough location of Tc. To do so we calculate (/p), the 
average of the thermal line (an unblocked Polyakov loop that winds once around the 
Euchdean time-torus), the lightest mass, arrit, that contributes to spatially separated 
correlations of such fines, and {Q^), the average value of the fiuctuations of the topo- 
logical charge Q. It is clear from Table |T5| that there is indeed a phase transition 
close to Lt = 5: the thermal line develops a non-zero vacuum expectation value, and 
consequently the lightest mass becomes ~ (the energy of the vacuum). We also see 
a very striking suppression of the topological susceptibihty, a'^Xt = {Q^)/{L^Lt) across 
Lt = 5. One observes a similar but much less dramatic behaviour in SU(3) |2^, and 
this suggests that Xt might be an order parameter for the deconfining phase transition, 
at least as A^ — > oo, with appropriate consequences for the axial U(l) anomaly. This 
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is a topic we shall expand upon elsewhere. 

Having established what is 'high' T in our calculation, we show in Table ^ the 
(screening) masses one obtains from spatially separated correlators oi k = 1 and k = 2 
loops that wind once around the spatial torus. At high T we expect a oc T^ (since g^ 
is now dimensionless) . We see that while our flux loop masses grow faster than T they 
grow less fast than T^. This should be no surprise; rather one should be surprised by 
the precociously early onset of the linear high-T behaviour in the case of D=2+l. 

From the flux loop masses we calculate the string tension in two ways: assuming 
no string correction, i.e. setting c^ = in eqn(^, and assuming a bosonic string 
correction. One might expect that at high T one should use a string correction that is 
half-way between, because one has lost one of the transverse dimensions. In any case it 
is clear that the high-T ratio is consistent with Casimir scaling but probably not with 
the MQCD formula. Just as for D=2+l we see a dip in the masses at T just below 
Tc. By contrast, for this lattice calculation, the low-T calculation on the 10^ lattice is 
consistent with MQCD but not really with Casimir scaling. 

Once again one might try to use dimensional reduction to relate the Casimir scaling 
we observe at high T in D=3-|-l to our observation of it, earlier on in this paper, at 
low T in D=2+l. However any such argument must address the caveats created by 
the presence of extra adjoint scalars after the reduction. We also note that Casimir 
scaling - at least for fc-strings and at very high T - has been predicted in a recent 
model calculation p3 which speculates that at high T there is a plasma of adjoint 



magnetic pseudoparticles 'dual' to the gluon plasma. (Again one might try to relate 



27| this to low-T near-Casimir scaling in D=2+l via dimensional reduction.) 



5.6 Unstable strings 

Our calculations have so far focussed on stable fc-strings. In addition to these, there are 
also unstable strings, which are energetically unfavourable. An unstable string should 
appear as a nearly stable excited state in the fc-string spectrum. It is clearly interesting 
to find out how the string tensions of these strings depend on their representations and 
on A^. We note that this is precisely the question addressed by the calculations in. 



for example, [§, |T5|, |T6|. In this Section we will investigate closed strings with the 
quantum numbers of strings connecting k quarks in a given irreducible representation 
for k = 2,3. 

The tensions of strings connecting sources in a given irreducible representation can 
be extracted from correlation functions of operators carrying the quantum numbers 
of that representation. For k = 2 the irreducible representations with two quarks 
are the symmetric and the antisymmetric representations. At A; = 3 the irreducible 
representations with three quarks are the totally antisymmetric, the totally symmetric 
and the mixed symmetry representation, which enters twice the decomposition of the 
tensor product. The general procedure to construct the relevant operators for a given 
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representation and the explicit form of the operators corresponding to the cases we 
shall investigate in this Section are given in Appendix |A.1| . 

At each A/'-ality (and at finite A^) only the string with the smallest tension is stable. 
For any k, this string is expected to be the string connecting k sources in the totally 
antisymmetric representation. Hence the smallest mass in the antisymmetric channel 
is related to the string tension at the given A/'-ality. Our calculations fulfill this expec- 
tation: the smallest mass extracted with the variational procedure and the smallest 
mass in the antisymmetric channel always agree well within errors, both in D=2+l 
and D=3+l. 

At fixed length an unstable fc-string is more massive than the stable string of the 
same A/'-ality. This is likely to give problems when looking at the exponential decay 
in time of correlation functions: the signal will decay too rapidly to allow a reliable 
extraction of masses. Indeed this is what happens in our D=3+l calculations: given the 
precision of our numerical data, it proves to be impossible to extract reliable masses. To 
overcome this problem, we should get closer to the continuum limit or use anisotropic 
lattices. We leave such a study for the future. Another crucial point is the overlap 
between the operators and the interesting states, which if the operators are constructed 
using standard techniques (as we have done in the present study) can get as bad as 0.5 
for unstable strings in D=3+l. (For comparison, still in D=3+l, the overlap between 
the stringy state corresponding to the antisymmetric representation and our operators 
is typically around 0.6.) Hence, to address questions connected to unstable strings 
a better overlap is required. This requires in turn an improvement of the standard 
smearing techniques. This is another problem we will investigate in the future. 

While we can not deal at the moment with unstable strings in D=3+l, our numer- 
ical results in D=2+l allow us to address the question there. In fact, our D=2+l 
calculations do not suffer from the same drawbacks as the D=3-|-l ones: our D=2+l 
results are accurate enough to see a clear exponential decay of the correlation functions 
over several lattice spacings and in D=2+l the overlap between the operators and the 
unstable strings is not worse than 0.85. Our numerical results for SU(4) and SU(6) are 



reported respectively in Tables [T^ and |18[ Apart from the symmetric representation 
of SU(6), for which we have masses for just two values of the lattice spacing, we can 
extrapolate the string tensions extracted from the masses listed in the Tables to the 
continuum limit by applying the same procedure used for stable strings. We then find 

o^^ I 2. 14 ±0.03 SU(4) 
i™ a \ 2.19 ±0.02 SU(6) ^^ ' 

for the string tensions in the k = 2 symmetric channels and 

lim ^^^^ = 2.71 ± 0.09 SU(6) (35) 

for the k = 3 mixed symmetry channel. (All these string tensions have been extracted 
using the bosonic string correction.) For the continuum value of the k = 3S string 
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tension, an estimate based on our data gives 

lim^^^^^ 3.72 ±0.12 SU(6). (36) 

Our numerical results can be compared with the predictions coming from Casimir 
scaling 

(Tk=2S cs I 2.40 SU(4) 

(38) 
(39) 



a \ 2.28... 


SU(6 


(7k=3M CS ^ ^ry 

a 


SU(6), 


^^^=^^^^3.85... 


SU(6). 



a 

We see that, at least for SU(6), these ratios satisfy approximate Casimir scaling, just 
like the stable /c-strings. As far as comparison with MQCD is concerned, we are not 
aware of calculations in that framework aimed to determine the string tensions of 
unstable strings. 

6 Discussion 

Our calculations in this paper were in two parts. In the first part we investigated 
directly the stringy nature of long flux tubes by calculating how the mass of a flux 
tube varies with its length /, and attempting to identify the 0(1//) term that is the 
leading string correction at large /. The coefficient of this term, cl, is directly related 
to the central charge of the effective string theory describing the long-distance physics 
of the confining fiux tube, and thus characterises its universality class. By considering 
flux tubes that wind around a spatial torus we were able to avoid the presence of 
explicit sources and the accompanying Coulomb term which can be so easily confused 
with the string correction. Working in SU(2) and at a reasonably small value of the 
lattice spacing a, we obtained in 3+1 dimensions a value cl = 0.98 ± 0.04 which is 
consistent with the simple bosonic string, for which cl = vr/3. In 2+1 dimensions we 
obtained cl = 0.558 ± 0.019, which is again consistent with the bosonic string value, 
which is Cl = it/6 in this case. In both dimensions our results would appear to exclude 
other plausible possibilities with, for example, some massless fermionic modes along 
the string. In addition, in the case of D=2+l our results were accurate enough to 
constrain the power of 1// to be unity (assuming it to be an integer) as one expects for 
an effective string theory. These results considerably increase, we believe, the evidence 
for the simple bosonic string model. There is however much scope for improving these 
calculations; not only in their accuracy and in the range of flux tube lengths studied, 
but also in exploring other values of a, so as to be confldent of the continuum physics, 
and in extending the calculations to other SU(iV) groups. 
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The second part of the paper dealt with /c-strings in SU(A^ > 4) gauge theories and, 
in particular, with the ratios of their string tensions, (Jk/u. Here we performed a range 
of calculations so as to be able to extrapolate to the continuum limit. In our D=3+l 
SU(4) and SU(5) calculations we found that the k = 2 string tension is much less than 
twice the fundamental k = 1 tension: the fc-string is 'strongly bound'. Moreover the 
values are consistent, at the 2a level, with both the M(-theory)QCD conjecture and 
with Casimir scaling (the two being numerically quite similar). In our SU(4) and SU(6) 
calculations in D=2+l we also found strongly bound /c-strings. However, although the 
calculated string tension ratios were again numerically close to both Casimir scaling 
and the MQCD formula, the results were accurate enough for us to see that the former 
works much better, and to observe deviations from both formulae. In addition to 
these continuum calculations we performed some finite temperature calculations at 
fixed a which showed that at high T, above the deconfinement transition, 'spatial' k- 
string tensions are consistent with Casimir scaling in both 2+1 and 3+1 dimensions. 
Moreover we found fairly convincing evidence, in D=2+l, for the approximate Casimir 
scaling of unstable strings. While it might be elegant if (approximate) Casimir scaling 
were to hold in D=3+l as well as in D=2+l, and at high T as well as at low T, the fact 
is that 3+1 dimensions may well differ from 2+1 dimensions, and it is important to 
perform calculations that are accurate enough to resolve between MQCD and Casimir 
scaling in D=3+l. Essentially this would require reducing our errors by a factor of 
two, an entirely feasible goal. 

We observed that near-Casimir scaling will arise naturally if the chromo-electric flux 
is homogeneous and the cross-section of the flux tube is (nearly) independent of the 
flux carried. We pointed out that the latter is not as implausible as it might at first 
seem: indeed it is what occurs in the deep-London limit of a superconductor. To ad- 
dress this possibility we performed some explicit numerical calculations of the fc-string 
width and these indicated that the width is indeed largely independent of k. The small 
k dependence that we did observe can, in principle, be related to the parameters of 
the dual superconductor, if such is the dynamics of confinement, and we intend to 
address this question elsewhere. There are a number of other interesting theoretical 
questions that this work suggests. How closely is the observed near-Casimir scaling 
of the ('spatial') fc-strings at high-T in D=3+l and at low T in D=2+l related by 
dimensional reduction? Equally, is the near-Casimir scaling at high-T in D=2+l, a 
simple reflection of the Casimir scaling of the linearly confining Coulomb interaction 
in D=l+1? This requires understanding whether the adjoint scalars, present after di- 
mensional reduction, significantly affect the string tension ratios. Another interesting 
question is how the string tension ratios, whether given by MQCD or Casimir scal- 
ing, reflect themselves in fc-vortex condensates, in the dual disorder loop approach to 



confinement [^], and whether this exposes any simple duality between Wilson loops 
and 't Hooft disorder loops. A calculation, illustrating how one might proceed, was 
outlined in P, |2^. A similar question can be posed in monopole models of confine- 
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ment, following upon the simple model calculations of higher charged string tensions 
29| , after Abelian projection, and in |^ with adjoint monopoles (at high T). A 



m 



quite different question is what are the implications of these tightly-bound /i;-strings 
for the mass spectrum of SU(A^) gauge theories. A simple and attractive model sees 
the glueball spectrum as arising from excitations of closed loops of fundamental flux 



30| . In such a model a non-trivial fc-string would provide a new sector of states whose 
masses are scaled up by a simple factor of o"jt/cr [^. The observation of something like 
this, when comparing the SU(3) and SU(4) spectra for example, would provide striking 
information on glueball structure. While our D=3+l mass spectrum calculations [Q 
are too crude to usefully explore this question, this is not the case in D=2-|-l (see e.g. 
m) and work on this question is proceeding. 



Note added: As this paper was being completed, a paper appeared [^] containing a 



calculation of A; = 1,2,3 string tensions in D=3+l SU(6) and addressing some of the 
questions addressed in this work. 
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Appendix 

This Appendix collects detailed proofs of some statements contained in the main ex- 
position. In Section [A.l| we will derive the explicit form of the operators carrying the 



quantum numbers of /c-strings. In Section [A.2| we will show how sources in a given rep- 
resentation can be screened by gluons. Finally, Section |A.3| will deal with the quadratic 



Casimirs of irreducible representations of SU(A^) and their relationship with Casimir 
scaling. Our calculations will be heavily based on Group Theory. In order to make 
this Appendix self-contained, we will recall some general results of Group Theory. For 
a wider introduction to the group theoretical background we refer to |3^ . 
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A.l Irreducible representations of SU(A/^) and k-stings 

SU(iV) is the group oi NxN unitary matrices. By definition, an object g* (i = 1, . . . , A^) 
transforms under tlie fundamental representation of SU(A^) if under the action of the 
group 

q'^^-^^Uiq\ (40) 

U G SU(A^) being the matrix that implements the transformation. 

The conjugated representation is related to the fundamental one by complex conjuga- 
tion. Following the standard notation, we indicate by g^ an object transforming under 
the conjugated representation. For A^ > 3 the fundamental and the conjugated repre- 
sentations are independent. In the following, we will call quarks objects transforming 
under the fundamental representation and antiquarks objects transforming under the 
conjugated representation of SU(A^). This terminology reflects the physics of QCD. 

Objects transforming under higher representations can be constructed from the ten- 
sor product of quarks and antiquarks, and their transformation laws can be easily 
deduced from the transformation law of the fundamental constituents. For instance 
{q ® qY^ = q^q^ under the action of SU(A^) transforms as follows: 

qW'''-^''^UlU^q'qK (41) 

The A/'-ality of a representation is defined as the number of quarks minus the num- 
ber of antiquarks modulo A^. AA-alities k < N/2 and N — k are related by complex 
conjugation. That operation corresponds to charge conjugation. 

The concept of A/'-ality at zero temperature is related to the symmetry under the 
centre of the gauge group, Z^: under such symmetry an object of A/'-ality k pick up a 
phase e^^'^*'^")/^, n = 0,1, . . . , N — 1. Since the centre symmetry at zero temperature 
is a good symmetry of the gauge theory and the gluons carry zero A/'-ahty, states with 
different A/'-ality cannot mix. 

In this paper, we are interested in the tensions of strings connecting sources with 
A/'-ality k > 1. Because of charge conjugation, the string tension associated to states 
of A/'-ality k < N/2 and N — k is the same. Hence we restrict ourself to fc < N/2, that 
is to say to states constructed from the tensor product of k quark^. At a given N the 
independent number of stable fc-strings is given by the integer part of N/2. 

As for the fundamental string, the string tension of a fc-string can be extracted by 
looking at the exponential decay of correlators of loop operators with the quantum 
numbers of that string. In order to identify the relevant operators, it is useful to 
decompose the tensor product of k quarks into irreducible representations^. To this 
purpose, the Young tableau technique can be used. 

^Here we are neglecting states with more than N quarks; this is correct as far as we are not 
interested in unstable strings or N is large enough. 

^In the following, even if we will omit for simplicity the word irreducible from time to time, we will 
consider only irreducible representations of the gauge group. 
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A Young diagram is a two-dimensional ensemble of boxes joined by one edge that 
respects the following rules: 

1. counting the rows from the top to the bottom, the number of boxes in the row i 
is greater than or equal to the number of boxes in the row j ii i < j; 

2. counting the columns from the left to the right, the number of boxes in the 
column i is greater than or equal to the number of boxes in the column j Hi < j. 

A valid Young tableau is for instance the following: 



In the Young tableau language, a quark is a single box, an antiquark is a column of 
N — 1 boxes and an object transforming under the adjoint representation [gluon) has 
A^ — 1 boxes in the first column and 2 boxes in the first row: 



N-1 



N-1 



quark 



antiquark 



gluon 



There is a one-to-one correspondence between Young diagrams and irreducible repre- 
sentations of SU(A^). Given a Young diagram with k boxes, the object transforming un- 
der the corresponding irreducible representation is constructed from the tensor product 
of k quarks by assigning an index to each box, symmetrising the product with respect 
to the indices that are on a given row for all rows and then antisymmetrising the result 
with respect to the indices that are on a given column for all columns. (Obviously, 
after the antisymmetrisation the result is no longer symmetric under permutation of 
indices on the same row.) 

The tensor product of two objects transforming under two given representations of 
the gauge group is constructed from the corresponding Young diagrams, A and B, 
according to the following rules: 

1. write down the two tableaux A and B labelling each box in the row i oi B hy i\ 
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2. starting from the first row of B, add the boxes of -B to A one-by-one in all the 
possible positions respecting the following rules: 

(a) the augmented diagram A' at each stage must be a legal Young diagram; 

(b) boxes with the same label must not appear in the same column of A'; 

(c) If we define at any given box position J numbers ni, . . . ,nj (J being the 
number of rows in B), each of them counting how many times the corre- 
sponding label of the boxes in B appears above and to the right of such a 
box, we must have rii > n2 > ■ ■ ■ > nj (this is to take into account the 
original antisymmetries of B); 

3. two diagrams with the same shape and the same labels are the same diagram; 

4. columns with A^ boxes must be canceled, since they correspond to the trivial 
representation of SU(A^). 

According to the above rules, the tensor product of two quarks decomposes as 



i.e. the irreducible representations of a state with two quarks are the symmetric and 
antisymmetric representations. For three quarks we have 



© 



where in addition to the symmetric and antisymmetric representation there is a repre- 
sentation with mixed symmetry entering twice the decomposition. 

The above results are the generalisation in SU(A^) of the familiar decompositions in 
SU(3) 3®3 = 6©3 and 30303 = 10 ©8©8©1. 

Once the symmetry of the states transforming under an irreducible representation has 
been worked out, it is easy to construct the operators implementing the transformation 
on such states, since those operators must have the same symmetry as the states on 
which they act. For the matrix elements oi k = 2 operators associated to strings 
connecting sources with two quarks we obtain 



^«m ~ o ( ^/ ^m 



urn: 



I ^ m I 1 



siL = Uuiui + ulK 



(42) 
(43) 
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while for k = 3 strings connecting three quarks we have 



A 



ijk 
Imn 



nijk 
Imn 



Imn 



UlUrnPn 



UlUtUi - UiUlJJ^ + UiutK + UfUlUi - UtUiK) , (44) 
i {uiUiU!; + UiUtUi + UlUlU'^ + UlutUl, + U^UlUi + U^Ui^Ui) , (45) 
\ (uiU^U!: - UtU^^ + UiUlUl - UiutUi) , (46) 



A, S and M being respectively the tensors corresponding to the antisymmetric, the 
symmetric and one of the two mixed symmetry representations. (The other mixed 
symmetry representation has k and j interchanged in eqn(^6|).) 
Taking the trace, we get 



TrA 

TrS 



I {{TrUy 



TrU' 



2 

- [{TrUy + TrU^ 



(47) 
(48) 



for k=2 and 



^ ({Trf/}3 - 3TrU{TrUy + 2TrU'^) , 


(49) 


^ ({Trf/}3 + 3TrU{TrUY + 2TrU^) , 


(50) 


\ {{TrUr - TrU') , 


(51) 



TrA 

TrS 
TrM 



for k=3. (The two different M's have the same trace.) 

After identifying U with the path ordered product of links around a non-contractible 
loop c that winds once around the spatial torus, Pc, we get that the relevant operators 
for k = 2 are TrP^ and {TrPc}^, while for A; = 3 we will be concerned with TrPj^, 
TrPc{TrPc}^ and {TrPc}^. These operators can be taken as a starting point for a 
variational procedure to extract the mass of flux tubes of A/'-ality k winding once 
around the periodic lattice, while studying directly the combination corresponding to 
a given irreducible representation is relevant in the context of unstable strings (see the 
following Section). 

The construction here explicitly provided for k = 2 and A; = 3 can be easily gener- 
alised to any k. 

A. 2 Gluon screening 

Since the gluons transform under a non-trivial representation of SU(A^), the interaction 
between them and the sources can change the original representation of the sources. 
The change of representation of the source is expected to renormalise the string tension 
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associated to the original representation and to make unstable heavier strings of given 
A/'-ality. 

Prom the point of view of Group Theory, the product of the interaction between 
sources and gluons transforms as the tensor product of the original representation and 
the adjoint representation. Consider for instance the following interaction: 



where the first tableau correspond to a source with k quarks and the second diagram 
is associated to a gluon. The interaction will produce an object transforming under a 
reducible representation. The irreducible representations entering the product of the 
interaction can be worked out according to the rules for the decomposition of a tensor 
product given in the previous section. Those representations fall into two categories: 
representations with k quarks and representations with N + k quarks. 

Let us consider the first case. In order to have a final state with k quarks, in the 
tensor product A^ boxes must be canceled (i.e. they are combined in such a way that 
they transform under the trivial representation). Given the diagram of a source with 
m boxes in the first column, the cancellation of N boxes in the tensor product requires 
that N — m boxes from the gluon are attached to the first column of the source. This 
can be done in two inequivalent ways: by taking the required objects all from rows 
other than the first one or by taking one box from the first row^. The possible ways of 
recombining the diagrams after the cancellation define possible representations of the 
interacting state. Those representations depend on the original representation of the 
sources, but not on A^. Similar considerations hold for the irreducible representations 
with N + k quarks entering the decomposition. 

An interacting state will be energetically favourable whenever it has a smaller string 
tension, so we expect that the interaction tends to transform the sources in a given 
representation of A/'-ality k to sources in the representation with the smallest string 
tension (which is the antisymmetric representation with k quarks in both the Casimir 
scaling and the MQCD scenario), i.e. that the gluons screen the sources down to the 
states with the smallest string tensionfl. However unstable strings are expected to be 



■^This argument should be refined if we were interested to the multiplicity with which each irre- 
ducible representation enters the decomposition of the tensor product. 

^Note however that there is a phase-space suppression factor for the screening of sources with a 
given number of quarks down to sources with a minor number of quarks inside the same A/'-ality class. 
Consider for instance the screening of a representation with N + k quarks to a representation with k 
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visible, since they should appear as nearly-stable excited states in the mass spectrum 
of the strings. 

Note that not all states with A/'-ality k are accessible to a given state. For instance, 
the interaction with one gluon does not allow to pass from the symmetric to the an- 
tisymmetric representation at A; = 3. However any state can be accessed by multiple 
interaction. 

A. 3 Quadratic Casimir operator and Casimir Scaling 

The quadratic Casimir operator of a representation TZ is defined as 

ct^^Y.T'^T'^ (52) 

a 

where the sum ranges over all the generators of the group in the given representation. It 
can be easily seen that the quadratic Casimir operator commutes with all the generators 
of the group. Hence, by virtue of the Schur's lemma, on a given representation it 
is proportional to the identity, i.e. it is identified by a number depending on the 
representation. We call such a number quadratic Casimir and we indicate it by C-n. If 
we normalise the trace of the identity in each representation to 1, we can write 

Cn = TrT''T\ (53) 

Cti can be easily computed starting from the Young tableau associated to the rep- 
resentation TZ as follows. For SU(A^), defines the A^-dimensional vectors 

Li = ^(iV- 1,-1,-1,-1,...,-!), 
L2 = ^(iV-2,iV-2,-2,-2,...,-2), 
U = ^(iV-3,iV-3,iV-3,-3,...,-3), 

Ln-1 = ^(1,1,1, ,l,...,-(iV-l)), 

2R = (iV-l,A^-3,A^-5,...,-(iV-3),-(iV-l)). 

With the vectors Lj, define 

N~l 
L=J2 ^^L^ (54) 



i=l 



quarks. Since the dimension of the former is proportional to N^^'^ ^ while the dimension of the latter 
is proportional to N'^ , there is a suppression due to the lack of final states proportional to 1/N^. 
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where Wi is given by the difference between the number of boxes in the row i and the 
number of boxes in the row i + 1 of the Young tableau. The quadratic Casimir is then 
given by 

Cn = \{L-L + 2R-L). (55) 

It is now easy to see that for an irreducible representation composed by k quarks the 
quadratic Casimir is given by the formula 

1 / m jL2\ 

Cn = i^ [Nk + ^n,(n, + 1 - 2^ - - j (56) 

where i ranges over the rows of the Young tableau (with m number of rows) and Ui 
is the number of boxes in the z-th row. For the antisymmetric and the symmetric 
representations of A/'-ality k we have 

c^^c,^-^^ ,5r, 

and 

^-^/^^. (5«) 

N"^ - 1 
Cf = (59) 

being the quadratic Casimir of the fundamental representation. 

For A; = 3 in addition to the symmetric and antisymmetric representations, there is 
(among others) the mixed symmetry representation, whose quadratic Casimir is 

Casimir scaling is the hypothesis that the string tension for a given representation 
is proportional to the quadratic Casimir. Hence, according to this hypothesis, at 
TV-ality k the smallest string tension is associated to sources in the representation 
with the smallest quadratic Casimir. By using eqn(|56|), it can be easily seen that 
the representation having the smallest quadratic Casimir is the totally antisymmetric 
representation composed by k quarks. To show this, let us prove as a preliminary step 
that if we increase the antisymmetries of a diagram keeping constant the number of 
boxes the quadratic Casimir decreases. In fact the difference between the quadratic 
Casimir of a given representation and of the representation obtained by moving a box 
of the original Young diagram from the j-th row to the h-th row with h > j is 

AC = nj-nh + h-j -l>0. (61) 

It is now easy to prove the main statement: at given number of boxes k < N/2, the 
antisymmetric representation is obtained from any given representation by iterating 
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the above procedure, with a series of steps where at each stage the quadratic Casimir 
decreases. Thus at given k the antisymmetric representation with k quarks has the 
smallest quadratic Casimir and for this reason the smallest string tension within the 
class of the representations with A/'-ality k in the Casimir scaling hypothesis. This 
fact does old even if we consider states with N + k quarks: the difference between the 
quadratic Casimir of the most antisymmetric representation with N + k quarks (which 
is the smallest at that number of quarks) and the totally antisymmetric representation 
with k quarks is 

AC = N -k. (62) 

Not all the representations with k quarks have a larger quadratic Casimir than a 
given representation with N + k quarks. For instance, the difference between the 
quadratic Casimir of the most antisymmetric representation with N + k quarks and 
the totally symmetric representation with k quarks is 

AC = N- e, (63) 

which is negative if k"^ > N. However, at large enough A^ and at a given k such 
a difference is positive and increases as A^. That is to say, we expect that the only 
relevant states of A/'-ality k in the large N limit are those composed by k quarks. 

The prediction of Casimir scaling for the ratio between the string tensions associated 
to states composed hj N + k quarks and the string tension of the fundamental repre- 
sentation in the limit A^ ^ oo is A; + 2. This result can be easily understood in terms 
of string counting: a state with N + k quarks can be seen as a possible state among 
those originated by the interaction between a state with k quarks and a gluon. The 
above result then tells us that at large A^ the energy of the composite state is equal to 
the sum of the energies of the constituents. 
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ami ; SU(2) ; D=2+l 



lattice 



8^64 
lOMS 
12236 
14236 
16^32 
20^32 
24232 
28^32 
32^32 
36^32 
40^32 



MC sweeps 



4x 10^ 
8 X 10^ 

10^ 
2 X 10*^ 
2 X 10*^ 
2 X 10*^ 
2 X 10^ 
2 X 10^ 
2 X 10^ 
2 X 10*^ 
2 X 10*^ 



p= 



0.1703(4) 
0.2167(5) 
0.2696(5) 
0.3219(6) 
0.3812(4) 
0.4917(7) 
0.5998(13) 
0.7101(10) 
0.8131(23) 
0.9175(33) 
1.0238(53) 



'safe' p = 



0.1703(4) 
0.2167(5) 
0.2696(5) 
0.3210(8) 
0.3806(6) 
0.4906(9) 
0.5998(13) 
0.7089(18) 
0.8131(23) 
0.9167(38) 
1.0238(53) 



low p 



0.1703(4 
0.2167(5 
0.2696(5 
0.3219(6 
0.3812(4 
0.4917(5 
0.6004(8 
0.7083(9 
0.8150(15) 
0.9195(19) 
1.0284(24) 



Table 1: The lightest mass ami of a fundamental string wrapped around a spatial 
torus. The first column comes from good fits to p = correlators, chosen so as to 
minimise the errors. The second column contains cautious 'very safe' estimates with 
larger errors. The third column uses both p = and, where useful, p ^ correlators. 



ct^ni,n;i 


3U(2) ; D=2+l 


L 


L' 


p= 


'safe' p = 


low p 


8 


10 


0.130(24) 


0.130(24) 


0.130(24) 


10 


12 


0.498(41) 


0.498(41) 


0.498(41) 


12 


14 


0.546(62) 


0.479(73) 


0.546(62) 


14 


16 


1.329(79) 


1.37(11) 


1.329(79) 


16 


20 


1.046(58) 


1.01(8) 


1.032(48) 


20 


24 


0.99(16) 


1.15(18) 


1.08(10) 


24 


28 


1.53(27) 


1.35(35) 


1.16(19) 


28 


32 


0.31(51) 


0.6(6) 


1.10(36) 


32 


36 


0.7(1.1) 


0.5(1.2) 


0.68(66) 


36 


40 


1.4(2.1) 


1.7(2.2) 


2.2(1.0) 



Table 2: Estimating the effective string correction coefficient from the masses of pairs 
of flux loops of length / = aL and /' = aL' respectively, using eqn(|l^). 
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Cs{l 


> lo) ; SU(2) 


D=2+l 






Lo 


p = 


xVdof 


'safe' p = 


xVdof 


low p 


xVdof 


14 


- 


- 


1.118(28) 


1.3 


1.104(18) 


1.6 


16 


1.071(26) 


1.2 


1.066(36) 


0.4 


1.070(22) 


0.4 


20 


1.091(63) 


1.4 


1.11(8) 


0.5 


1.105(42) 


0.3 


24 


1.02(16) 


1.9 


1.01(16) 


0.5 


1.12(9) 


0.4 


28 


0.52(28) 


0.2 


0.69(35) 


0.1 


1.10(17) 


0.6 


32 


0.9(7) 


0.1 


0.88(75) 


0.2 


1.21(40) 


1.0 



Table 3: Estimating the string correction coefficient from a fit of eqn(^ to the masses 
of all the flux loops of length I > Iq = oLq. In each case we show the x^ /dof of the 
best fit. 



ami ; SU(2) ; D=3+l 


lattice 


MC sweeps 


p = 


'safe' p = 


low p 


'safe' low p 


10^60 


10^ 


0.1679(14) 


0.1679(14) 


0.1679(14) 


0.1679(14) 


12^48 


2x 10^ 


0.2073(14) 


0.2073(14) 


0.2073(14) 


0.2073(14) 


14^36 


4x 10^ 


0.2632(13) 


0.2606(16) 


0.2636(12) 


0.2606(16) 


16^28 


6 X 10^ 


0.3230(18) 


0.3230(18) 


0.3302(11) 


0.3230(18) 


20^ 


6x 10^ 


0.4468(15) 


0.4416(23) 


0.4408(20) 


0.4408(20) 


24^ 


8x 10^ 


0.5476(22) 


0.5476(22) 


0.5459(15) 


0.5459(15) 


32^24 


4x 10^ 


0.7598(75) 


0.7469(115) 


0.7549(50) 


0.7496(58) 



Table 4: The lightest mass ami of a fundamental string wrapped around a spatial 
torus. The ffist column comes from good fits to p = correlators, chosen so as to 
minimise the errors. The second column contains cautious 'very safe' estimates with 
larger errors. The third and fourth columns use both p = and, where useful, p 7^ 
correlators. 





c^^^(/,/');SU(2); 


D=3+l 




L 


V 


p = 


'safe' p = 


low p 


'safe' low p 


10 


12 


0.15(6) 


0.15(6) 


0.15(6) 


0.15(6) 


12 


14 


0.79(8) 


0.69(9) 


0.81(8) 


0.69(9) 


14 


16 


1.11(12) 


1.26(13) 


1.44(9) 


1.26(13) 


16 


20 


1.46(9) 


1.28(11) 


0.95(8) 


1.26(10) 


20 


24 


0.60(15) 


0.92(18) 


0.88(15) 


0.88(15) 


24 


32 


1.17(32) 


0.66(47) 


1.06(21) 


0.85(24) 



Table 5: Estimating the effective string correction coefficient from the masses of pairs 
of flux loops of length / = aL and /' = aV respectively, using eqn(p!l|). 
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Cs{l > /o) ; SU(2) ; D=3+l 


Lo 


'safe' p = 


xVdof 


low p x^/dof 


'safe' low p 


x'/dof 


14 
16 
20 


1.19(4) 

1.14(6) 

0.88(16) 


1.5 
1.7 
0.2 


0.94(4) 0.2 
0.95(12) 0.4 


1.150(32) 
1.087(52) 
0.87(11) 


2.2 
2.2 
0.0 



Table 6: Estimating the string correction coefficient from a fit of eqn(^ to the masses 
of all the flux loops of length I > Iq = aL^. In each case we show the x^ /dof of the 
best ffi. 



D=3+l 


SU(4) 


su 


.5) 


13 


a^a 


/3 


a^a 


10.55 


0.372 


16.755 


0.384 


10.70 


0.306 


16.975 


0.303 


10.90 


0.243 


17.270 


0.245 


11.10 


0.202 


17.450 


0.222 


11.30 


0.170 


- 


- 



Table 7: Setting the scale of a: the string tension for our SU(4) and SU(5) lattice 
calculations in D=3+l. 



D=3+l ; SU(4) 


/? 


lattice 


MC sweeps 


amfc=i 


amk=2 


10.55 


8^ 


2x 10^ 


0.973(17) 


1.456(30) 


10.70 


6^16 


5x 10^ 


0.268(8) 


0.329(12) 


10.70 


8^12 


10^ 


0.564(10) 


0.763(24) 


10.70 


10^ 


10^ 


0.8375(92) 


1.197(18) 


10.70 


12^ 


10^ 


1.033(11) 


1.456(37) 


10.70 


144 


10^ 


1.201(34) 


1.780(60) 


10.70 


16^ 


10^ 


1.316(78) 


2.35(27) 


10.90 


12^ 


10^ 


0.622(7) 


0.896(11) 


11.10 


16^ 


10^ 


0.585(8) 


0.836(21) 


11.30 


20^ 


10^ 


0.5265(56) 


0.740(16) 



Table 8: Masses of the k = 1 and k = 2 flux loops that wind once around the spatial 
torus. In D=3+l SU(4), for the lattices and couplings shown. 
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D=3+l ; SU(5) 


/9 


lattice 


MC sweeps 


amk=i 


amk=2 


16.755 
16.975 
17.270 
17.450 


8^ 
10^ 
12^ 
16^ 


10^ 
2.0 X 10^ 
1.4 X 10^ 

10^ 


1.051(13) 
0.816(12) 
0.638(9) 
0.723(10) 


1.70(7) 
1.239(51) 
1.061(28) 
1.168(62) 



Table 9: Masses of the k = 1 and k = 2 flux loops that wind once around the spatial 
torus. In D=3+l SU(5), for the lattices and couplings shown. 



D= 


2+1 


SU(4) 


SU(6) 


P 


a^/a 


P 


a^/cr 


18.0 


0.442 


42.0 


0.436 


21.0 


0.361 


49.0 


0.353 


28.0 


0.252 


60.0 


0.274 


33.0 


0.208 


75.0 


0.211 


45.0 


0.147 


108.0 


0.141 


60.0 


0.108 


- 


- 



Table 10: Setting the scale of a: the string tension for our SU(4) and SU(6) lattice 
calculations in D=2+l. 







D=2+l ; 


SU(4) 


p 


lattice 


MC sweeps 


amk=i 


amk=2 


18.0 


83 


2x 10^ 


1.497(18) 


2.022(66) 


21.0 


103 


2x 10^ 


1.223(11) 


1.690(38) 


28.0 


4^48 


10^ 


0.1347(12) 


0.1836(21) 


28.0 


6^28 


10^ 


0.2724(15) 


0.3630(26) 


28.0 


8^16 


10^ 


0.4276(33) 


0.5652(64) 


28.0 


10^16 


1.5 X 10^ 


0.5720(30) 


0.7824(63) 


28.0 


12^ 


2x 10^ 


0.7152(39) 


0.9718(95) 


28.0 


163 


2x 10^ 


0.9937(58) 


1.345(14) 


33.0 


163 


2x 10^ 


0.6622(29) 


0.914(6) 


45.0 


243 


2x 10^ 


0.4974(22) 


0.6815(40) 


60.0 


323 


2x 10^ 


0.3571(14) 


0.4888(15) 



Table 11: Masses of the k = 1 and k = 2 flux loops that wind once around the spatial 
torus. In D=2+l SU(4), for the lattices and couphngs shown. 
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D= 


2+1 ; SU(6) 






/9 


lattice 


MC sweeps 


amk=i 


amk=2 


amk=3 


42.0 


8^ 


3x 10^ 


1.453(14) 


2.51(10) 


2.85(30) 


49.0 


10^ 


2x 10^ 


1.194(9) 


2.030(42) 


2.21(9) 


60.0 


12^ 


2x 10^ 


0.8575(47) 


1.443(15) 


1.621(25) 


75.0 


16^ 


2x 10^ 


0.6825(31) 


1.1305(72) 


1.275(11) 


108.0 


243 


3x 10^ 


0.4552(13) 


0.7534(21) 


0.8424(98) 



Table 12: Masses of the k = 1, k = 2 and A; = 3 flux loops that wind once around the 
spatial torus. In D=2+l SU(6), for the lattices and couplings shown. 





D=2+l ; SU(4) 


; /9 = 28 




amk=i{L) 


amk=2{L) 


L± 


L = 8 


L = 12 


L = 8 


L = 12 


2 


1.188(10) 


1.811(23) 


1.552(28) 


2.27(10) 


3 


0.738(9) 


1.1900(73) 


1.008(11) 


1.45(15) 


4 


0.5093(37) 


0.819(13) 


0.667(9) 


1.04(4) 


5 


0.4297(36) 


0.6814(83) 


0.5583(62) 


0.899(12) 


6 


0.4207(28) 


0.6965(46) 


0.5504(57) 


0.869(22) 


8 


0.4276(33) 


0.710(5) 


0.5652(64) 


0.947(9) 


10 


0.4274(30) 


0.706(5) 


0.5770(56) 


0.971(11) 


12 


0.4327(30) 


0.715(4) 


0.5985(58) 


0.972(10) 


16 


0.4347(28) 


0.721(5) 


0.598(7) 


0.948(24) 


20 


0.4331(30) 


0.713(4) 


0.600(5) 


0.991(12) 



Table 13: Masses of flux loops that wind once around a spatial torus of length L, as 
a function of the length, Lj_, of the other spatial torus. 
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D=2+l ; SU(6) ; /5 = 60 




amk=i{L) 


amk=2{L) 


amk=3{L) 


l± 


L = 10 


L = 12 


L = 10 


L = 12 


L = 10 


L = 12 


2 


1.40(13) 


2.055(8) 


2.66(20) 


2.70(41) 


3.11(4) 


2.63(91) 


3 


1.050(29) 


1.227(56) 


1.785(38) 


2.069(51) 


1.929(70) 


2.33(21) 


4 


0.692(12) 


0.841(18) 


1.178(14) 


1.419(28) 


1.360(27) 


1.611(54) 


5 


0.6767(66) 


0.8455(65) 


1.064(14) 


1.280(11) 


1.163(16) 


1.389(31) 


6 


0.6915(67) 


0.8462(59) 


1.084(12) 


1.368(23) 


1.200(21) 


1.531(31) 


8 


0.7005(52) 


0.8671(54) 


1.137(14) 


1.390(23) 


1.194(65) 


1.590(37) 


10 


0.7015(49) 


0.843(13) 


1.180(14) 


1.435(16) 


1.327(27) 


1.658(55) 


12 


- 


0.8575(47) 


- 


1.443(15) 


- 


1.621(25) 



Table 14: Masses of flux loops that wind once around a spatial torus of length L, as 
a function of the length, L_i_, of the other spatial torus. 



D=3+l ; SU(4) ; finite T 


Lt 


T/> 


(Ip) 


aiTLt 


m 


2 
3 

4 
5 
6 


1.63 
1.09 
0.82 
0.65 
0.54 


0.5669(1) 
0.3414(1) 
0.1922(2) 
0.0041(40) 
-0.0006(10) 


0.000 

0.003 

0.005 

0.052(1) 

0.230(8) 


0.002(1) 
0.150(15) 
1.363(53) 


10 


'0.33' 


- 


0.838(9) 


2.56(9) 



Table 15: Calculations at finite temperature, T, on 10'^ x Lt lattices. In SU(4) at 
(3 = 10.7. We show the thermal line average, (Ip), the lightest mass coupled to the 
thermal line, arrit, and the size of the topological fluctuations, (Q^). 
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D=3+l ; SU(4) ; finite T 




0"fe=2/0- 


Lt 


amk=i 


amk=2 


Cs = 


Cs = l 


2 
3 

4 
5 
6 


2.71(17) 
1.523(28) 
1.044(9) 
0.800(9) 
0.727(16) 


2.7(7) 
2.11(10) 

1.40(3) 
1.101(20) 

0.90(4) 


1.385(70) 
1.341(30) 
1.376(30) 
1.238(55) 


1.361(66) 
1.310(28) 
1.333(26) 
1.208(54) 


10 


0.838(9) 


1.197(18) 


1.429(27) 


1.382(23) 



Table 16: Calculations at the finite temperatures listed in Table ^. We list the 
'spatial' loop masses and the corresponding ratio of 'spatial' string tensions, calculated 
with and without a bosonic string correction. 



D=2+l ; SU(4) 


P 


lattice 


amk=2s 


28.0 
33.0 
45.0 
60.0 


123 
163 
243 
323 


1.601(23) 

1.352(15) 

0.979(40) 

0.7694(70) 



Table 17: Masses of fiux loops in the k = 2 symmetric representation. 
SU(4), for the lattices and couplings shown. 



In D=2+l 



P 



49.0 
60.0 
75.0 
108.0 



D=2+l ; SU(6) 



lattice 



103 
123 
163 
243 



amk=2s 



2.58(19) 
1.928(41) 
1.512(15) 
1.0210(47) 



amk=3M 



2.44(12) 

1.845(37) 

1.274(25) 



amk=3S 



2.270(87) 
1.646(17) 



Table 18: Masses of fiux loops in the k = 2 symmetric {k = 2S), k = 3 mixed 
[k = 3M) and k = 3 symmetric {k = 3S) representations. In D=2+l SU(6), for the 
lattices and couplings shown. 
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Figure 1: The D=2+l effective string correction coefficient estimated from tlie masses 
of ffux loops of different lengtlis (indicated by the span of the horizontal error bar) 
using eqn (|Tl|) . The solid line is what one expects for a simple bosonic string. 
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Figure 2: The D=2+l string correction coefficient estimated by fitting the masses of 
all flux loops with length greater than L, as a function of L. 
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Figure 3: The D=3+l effective string correction coefficient estimated from the masses 
of fiux loops of different lengths (indicated by the span of the horizontal error bar) 
using eqn(|l^). The solid line is what one expects for a simple bosonic string. For 
comparison the dashed line indicates the value for the Neveu- Schwartz string. We use 
masses from the third column in Table ^. 



49 



1.5 



X> 



o 

f 



0.5 



12 



16 



20 



Figure 4: The ratio oi k = 2 and k = 1 string tensions in D=3+l SU(4) at /3 = 10.7 
extracted from flux loops of length / = aL. We show values extracted using a bosonic 
string correction, (•), and no string correction at all (o). 
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Figure 5: The ratio of A; = 2 and k = \ string tensions in our D=3+l SU(4) (•) 
and SU(5) (o) lattice calculations plotted as a function of a^o. Extrapolations to the 
continuum limit, using a leading Oia^^ correction, are displayed. 



51 



1.5 

ffl m ^ ® ® ffl 

® fD -• • • 



^ 

d" 



0.5 



I . 1 . 1 . 1 . 1 . 

4 8 12 16 20 



Figure 6: The ratio of /c = 2 and k = 1 string tensions in D=2+l SU(4) at /9 = 28.0 
extracted from flux loops of length / = aL. We show values extracted using a bosonic 
string correction, (•), and no string correction at all (o). 
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Figure 7: The ratio of /c = 2 and k = 1 string tensions in our D=2+l SU(4) (•) and 
SU(6) (o) lattice calculations plotted as a function of a^o". Also shown is the k = 3 to 
k = 1 ratio (o) in SU(6). Extrapolations to the continuum limit, using a leading 0{a'^) 
correction, are displayed. 
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Figure 8: The masses of the k - 
D=2+l SU(4) gauge theory at [3 



1 and k = 2 flux loops of length L = 8 in the 
28, versus the size of the transverse spatial torus 



-'perp 



L±. Shown is the dependence in eqn(^3|) fitted to the smallest values of L 
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Figure 9: The ratio of the k = 2 to k = 1 string tensions in the D=2+l SU(4) gauge 
theory at /9 = 28 versus the size of the transverse spatial torus Lperp = L±. With (•) 
and without (o) a (bosonic) string correction. 
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Figure 10: The masses of the k = 1, k = 2 and /c = 3 flux loops of length L = 10 in the 
D=2+l SU(6) gauge theory at /5 = 60 versus the size of the transverse spatial torus 
Lperp = L±. Shown is the dependence in eqn(p3D fitted to the smallest values of Lperp- 
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Figure 11: The ratio of the k = 2 to k = 
in the D=2+l SU(6) gauge theory at P 



1 (•) and k = 3 to k = 1 (o) string tensions 
= 60 versus the size of the transverse spatial 



torus Lperp = L±. The (bosonic) string correction has been included. 
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